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Abstract. 

Electrostatic correlations play an important role in physics, chemistry and biology. 
In plasmas they result in thermodynamic instability similar to the liquid-gas phase 
transition of simple molecular fluids. For charged colloidal suspensions the electrostatic 
correlations are responsible for screening and colloidal charge renormalization. In 
aqueous solutions containing multivalent counterions they can lead to charge inversion 
and flocculation. In biological systems the correlations account for the organization 
of cytoskeleton and the compaction of genetic material. In spite of their ubiquity, the 
true importance of electrostatic correlations has become fully appreciated only quite 
recently. In this paper, I will review the thermodynamic consequences of electrostatic 
correlations in a variety of systems ranging from classical plasmas to molecular biology. 
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1. Introduction 

Although the liquid state theorists have become quite accustomed to look at the 
correlation functions for simple and complex fluids, many of the thermodynamic 
consequences of correlations have not been fully appreciated. To some extent this is the 
result of the success of mean-field theories for simple fluids. The strategy of applying the 
mean-field approximation to the many body problems goes all the way to the pioneering 
work of van der Waals on liquid-gas phase separation The success of this theory 
and its physical transparency has set the stage for future applications of mean-field 
ideas. These came in the form of the Curie- Weiss theory of magnetism [H and the 
Gouy-Chapman ^ theory of diffuse ionic layers. 

There are, however, some very familiar systems in which the mean-field contribution 
to the free energy is identically zero. One such system is the classical two component 
plasma of positive and negative ions. In order for the thermodynamic limit to exist, 
charge neutrality constraint must be imposed. However, for a bulk charge-neutral system 
the average electrostatic potential is zero, which means that the mean-field contribution 
to the total free energy vanishes. Thus, the electrostatic free energy of a two component 
plasma is entirely due to positional correlations between the positive and negative ions. 
At low temperatures these correlations become so strong as to lead to phase transition 
in which the plasma separates into two coexisting high and low density phases |^ . 

Polar fluids provide another example of a system in which the electrostatic 
correlations strongly affect the thermodynamics. Perhaps the simplest model of a polar 
fluid is a system of dipolar hard spheres (DHS). The phase structure of DHS is 
quite interesting and deserves a separate review [0, H, |^. Here we shall confine our 
attention to the low density disordered fluid phase. Since the average electric field 
inside a dipolar fluid is zero, it is evident that the mean-field contribution to the free 
energy also vanishes, and all the non-trivial thermodynamics is, once again, the result 
of electrostatic correlations. 

The thermodynamics of DHS is particularly tricky because of the unscreened long 
range interactions. In fact the very existence of the thermodynamic limit for DHS 
has been proven only recently Nevertheless, it has been taken for granted that 



if the temperature is sufficiently low, the DHS will phase separate into coexisting 
liquid and gas phases. Indeed, all the theories have been predicting exactly this kind of 



behavior |TT], |T^. It came, therefore, as a great surprise when the simulations in the early 
90's failed to locate the anticipated liquid-gas critical point |T3|, |T5[. Instead as the 
temperature was lowered, the simulations found chains of aligned dipoles. Formation of 
weakly interacting chains, a consequence of strong positional and directional correlations 
between the dipolar particles, prevented the liquid-gas phase separation from taking 



place [|TB|, [T^, ig |ra 

Electrostatic correlations are also crucial in charged colloidal suspensions |2^, [21 
In these systems the correlations come on two different levels. First, there are 
very strong positional correlations between the highly charged colloidal particles and 
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their counterions. These correlations lead to charge renormalization and to 
screening of Coulomb interactions between the colloidal particles. In water with 
monovalent counterions, charge renormalization stabilizes colloidal suspension against 



phase separation |^3|, In the presence of multivalent counterions, however, the 

layers of condensed counterions on different colloids can become strongly correlated, 
leading to a net attraction between the like-charged colloids [E^, EBI and to the phase 



separation ^7\, Pq] . 

A similar kind of behavior was also observed in a number of important biological 
systems. Thus, it was noted that like-charged macromolecules can attract each other 
in solutions containing multivalent counterions. This attraction manifests itself in in 
vitro formation of toroidal aggregates of concentrated DNA |2^, |3^, similar to the one 



found in bacteriophage heads and in bundle formation of F-actin and tobacco 
mosaic virus |3^. A number of models have been suggested to explain these curious 
phenomena. The fundamental ingredient in all of these models are the electrostatic 
correlations |3^, |3^, |35| . 

Although the phenomena described above are quite complex, we can get a long 
way towards understanding them by considering some surprisingly simple models and 
theories. In fact, we shall demonstrate that a lot of the physics of electrostatic 
correlations is contained within the Debye-Hiickel theory |^ introduced 80 years ago 
as a way of accounting for the unusual thermodynamic properties of strong electrolytes. 

Consideration of only simple physical theories in this review is partially pedagogical, 
designed for a broad audience not necessarily familiar with the complex machinery 
of correlation functions and field theories of the modern statistical mechanics. For 
Coulomb systems there is, however, an additional benefit. It is often found that the 
more sophisticated theories fail when applied to strongly correlated charged fluids. For 
example, the field theoretic calculations of Netz and Orland find that for charge 
asymmetric {z : 1) electrolytes the reduced critical temperature is a strongly increasing 
function of charge asymmetry. A completely opposite behavior is observed in computer 
simulations, the critical temperature decreases and the critical density increases with the 
charge asymmetry. In fact, the field theoretic predictions for the critical temperature 
of asymmetric electrolytes are so far-off that they have to be divided by a factor of 
six just to make them fit on the same graph (Fig. |ID with the results of simulations 
and of concurrent theories! The dramatic failure of field theoretic calculations can be 
attributed to their intrinsically perturbative nature. Similarly the integral equations, 
which have proven to be very successful for simple molecular fluids, fail to even converge 
for strongly asymmetric electrolytes. Furthermore, it has been known for a long time 
that the Hypernetted Chain Equation (HNC) which is often used to study the Coulomb 
systems [^, does not posses a true critical region, but only a "no solution zone" on 
the border of which compressibility goes to zero with a square root singularity |^ . 
This is a completely wrong behavior, since the compressibility must diverge at the 
critical point. All these should be contrasted with the physically based Debye-Hiickel- 
like theories, which are in qualitative and often in quantitative agreement with the 
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Figure 1. Estimates of tlie reduced critical temperature, T*, and density, p*, for 
the (z : 1), cliarge- asymmetric (but equisized) primitive model showing, as labeled, 
the predictions of pure Debye-Hiickel theory (without hard cores), (b) the MSA, (c) 
the SPB approximation (d) the Netz-Orland field-theoretic treatment which, 
for T* ^ have been divided by factors of 6 and 12 in order to bring them on to the 
plot (see labels NO' and NO" respectively), and (e) the DHBjCI theory (following the 
Fisher-Levin approach (5[ ^), point z = 3 is a preliminary result. The large open 
circles for z = 2 and 3 represent, the Monte Carlo simulations of Panagiotopoulos and 
Fisher while the estimate for T*(4) follows from Camp and Patey Q. (After S. 
Banerjee and M. E. Fisher, to be published). 



simulations and experiments, Fig. |1|. For some important many body systems, however, 
the integral equations provide the most accurate results. For example, predictions for 
the electrostatic free energy of the one component plasma obtained using the HNC 
equation are in excellent agreement with the Monte Carlo simulations The integral 
equations were also the first to account for the correlation induced attraction between 
the like-charged macromolecules I^S], In general, as long as one stays away from 

the phase transitions, the integral equations provide one of the sharpest tools available 
to a statistical physicist or chemist. Unfortunately, the approximations involved in 
constructing the integral equations are not very clear. There exists a great variety 
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of closures to the Ornstein-Zernike equation, each working well for specific kind of 
problems. Because of their complexity, I will not talk about the integral equations in 



this review, referring the interested reader to the literature [E? 



2. Electrolyte solutions 

Since the work of Faraday early in the 19'th century the flow of electricity has been 
associated with the movement of charged particles. The nature of ions (from Greek 
meaning wanderers), however, was not established. It was Arrhenius who in 1887, 
following the experimental work of van't Hoff on osmotic pressure of electrolyte solutions, 
proposed that when salts and acids are dissolved in water they become ionized ||48[] . 



Arrhenius suggested that NaCl dissociates forming cations Na^ and anions Cl~ . In the 
spirit of the mean-field theory introduced earlier by van der Waals [|l| , Arrhenius argued 
that since the anions and the cations are on average uniformly distributed throughout 
the solution, the average electrostatic potential inside the electrolyte is zero. He then 
concluded that on average there should not be any interaction between the ions, and 
the osmotic pressure of, say, IM solution of NaCl should be equivalent to the osmotic 
pressure of 2M solution of a non-electrolyte. All the deviations from this simple rule 
Arrhenius attributed to the incomplete dissociation of electrolyte, which he then treated 
as a problem of chemical equilibrium, the thermodynamics of which has already been 
developed by Gibbs ten years earlier |^9[. Soon, however, it became clear that while 



the theory was working well for weak electrolytes, such as Br0nsted acids and bases, it 
failed for strong electrolytes such as NaCl and HCl, which remained fully ionized even 
at fairly large concentrations. The disagreements between the theory and the experiment 
could not be accounted for by the postulate of incomplete dissociation. It appeared that 
there was a fundamental flaw in the theory advanced by Arrhenius, which relied on the 
mean-field assumption of non-interacting ions. The situation remained unclear for about 
30 years, with various proposals made on how to incorporate the ionic interactions into 
the framework of Arrhenius' theory None of these have proven very successful 



at explaining the experimental measurements, until Debye and Hiickel published their. 



now famous, theory of strong electrolytes [36|. The fundamental insight of Debye and 



Hiickel {DH) was to realize that although ions on average are randomly distributed, 
there exist strong positional correlations between the anions and cations. The depth of 
Debye's insight can be judged from the fact that he understood the role of electrostatic 
correlations significantly before the correlation functions became the standard tool of 
working physicists. Since so many of our results will be based on the fundamental ideas 
of Debye and Hiickel, their theory will provide the starting point for our discussion of 
thermodynamics of electrostatic correlations. 
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Figure 2. The configuration of closest approacfi between two oppositely charged ions. 
The dashed curve delimits the region into which no ions can penetrate, due to the hard 
core repulsion. 



2.1. The Debye-Hiickel theory 

Consider the simplest model of an electrolyte solution confined to volume V. The 
ions will be idealized as hard spheres of diameter a carrying charge ±g at their centers. 
The charge neutrality of solution requires that = N_ = N/2. The solvent will be 
modeled as a continuum of dielectric constant e. Although the average potential inside 
the electrolyte is zero, there are strong positional correlations between the oppositely 
charged ions. It is convenient to work in spherical coordinates. To calculate the 
correlational contribution to the Helmholtz free energy, let us fix one ion of charge 
+q at the origin r = and see how the other ions will distribute around it, see Fig. 
Inside the region < r < a there are no other charges except for the one fixed at the 
origin, and the electrostatic potential 0(r) satisfies the Laplace equation, 

V^<j) = 0. (1) 

For r > a the electrostatic potential satisfies the Poisson equation 

4-7r 

= -^P,(r) , (2) 

where the charge density can be expressed in terms of the charge-charge correlation 
functions g++{r) = g (r) and g+-{r) = g^+{r) 

Pq{r) = qp+g++{r) - qp-g+-{r) . (3) 

The average densities of positive and negative ions are p+ = N^/V, p_ = N_/V; 
P+ = P- = P/2- 

The correlation functions can be written in terms of the potential of mean force Wij 

grAr) = e-^--^'-) , (4) 

where (3 = l/ksT. The Wij{r) is the work required to bring ions i and j from infinity 
to separation r inside the electrolyte solution. In their paper Debye and Hiickel made 
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an implicit approximation of replacing the potential of mean force by the electrostatic 
potential 

Wij{r) = qj(f)i{r) , (5) 

where qj is the charge of j'th ion and </)j(r) is the electrostatic potential at distance r 
from the ion i fixed at the origin r = 0. With this approximation, Eq.(H) reduces to the 
non- linear Poisson-Boltzmann equation (PB), 

= [qp+e-^'^^ - gp_e+^'"^] = ^ sinh(/5g0) . (6) 

Debye and Hiickel proceeded to linearize this equation. Technically, linearization is only 
valid if l3q(f) <^ 1, however, being practically minded Debye and Hiickel linearized first 
and worried about the consequences later. As was noted later by Onsager, linearization 
of Eq.(|^) is a necessary step in order to produce a self-consistent theory |]SD[- The 
linearized Poisson-Boltzmann equation reduces to the Helmholtz equation 

= , (7) 
where the inverse Debye length is 



c-i - /47rg2p 

^-='^ = V^' 

The Laplace equation (|l|) for r < a and the Helmholtz equation (J^) for r > a 
must be integrated, subject to the boundary condition of continuity of the electrostatic 
potential and the electric field across the boundary surface r = a. For r < a the 
electrostatic potential is found to be 

<P<{r) = 7^ r , (9) 

er e(l + naj 

while for r > a, 

, , q9(Ka)e~'^'^ ^, e.^ , 

Mr) = , 0{x) = ——- . 10 

er (1 -f X) 

Equation ( [TI1| ) shows that the electrostatic potential produced by the central charge is 

exponentially screened by the surrounding ionic cloud. Because of the hardcore repulsion 

the screening, however, appears only at distances larger than r = a. This accounts 

for the presence of geometric factor 6{Ka) in Eq. ([T0|). The screening of electrostatic 

interactions inside the electrolyte solutions and plasmas is responsible for the existence 

of thermodynamic limit in these systems with extremely long range forces. 

The electrostatic potential </'<(r), Eq. (P), consists of two terms: the potential 

produced by the central ion q/er, and the electrostatic potential induced by the 

surrounding ionic cloud, 

e(l + Ka) 

The electrostatic free energy can now be obtained using the Debye charging process in 
which all the ions are simultaneously charged from zero to their full charge, 

F^^ = Nq [ dXil){\q) . (12) 

JQ 
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The calculation is very similar to the one used to obtain the electrostatic energy 
stored in a capacitor. While performing the charging it is important to remember that 
/t(Ag) = Xn^q). Defining the free energy density as / = F/V, the integral in Eq. (0) 
can be performed explicitly yielding 



i2 



\n{Ka + 1) — Ka 



(13) 



V 47ra3 

For large dilutions Eq. (|l^) reduces to the famous Debye limiting law, 

Given the free energy, the limiting laws for the osmotic pressure and activity can be 
easily found 0. 

The free energy is not analytic at p = 0. The singularity at p = is a consequence 
of long-range Coulomb interactions, which also manifest themselves in the divergence 



of the standard virial expansion ||51|. The total free energy of the electrolyte F is the 
sum of electrostatic Eq. (pTS]), and entropic contributions. The entropic contribution to 
the free energy arises from the integration over the momentum degrees of freedom in 
the partition function, and is equivalent to the free energy of an ideal gas, 

ppent ^ \n[p+A^] -N+ + N_ ln[p_A^] - = A^ln[pAV2] - N , (15) 
where the de Broglie thermal wavelength is 

A = ^ ^ . (16) 

The osmotic pressure of the electrolyte is 
dF 



^ dV 



dF 

p 



V op 



which can also be expressed in terms of the Legendre transform of the negative free 
energy density — / 0, 

P = -/ + pp, (18) 

where the chemical potential is 

_ 

dN 

It is a simple matter to see that below the critical temperature Tc the total free 
energy F = irent_|_^ez ^g^jjg ^ convex function of the electrolyte concentration. This 
implies the presence of a phase transition. Alternatively the phase separation can be 
observed from the appearance of a van der Waals loop in the osmotic pressure Eq. ([I8|) , 
below the critical temperature T^. The critical parameters are determined from 

dP , , 

g^=0. (20) 

1^ = 0- (21) 
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The coexistence curve can be obtained using the standard Maxwell construction. It 
is convenient to define the reduced temperature and density as T* = fc^Tae/g^ and 
p* = pa? . The critical point of the plasma, within the DH theory, is found to be located 
at [|, I 

= ^ , (22) 



16 



and 



It is interesting to note that at criticality k = 1/a. This means that in spite of a very low 
concentration of electrolyte at the critical point, the screening remains very strong. We 
also observe that the reduced critical temperature for the electrolyte is almost an order 
of magnitude lower, than for systems in which the particles interact by the short-ranged 
isotropic potentials. Since the critical point within the DH theory occurs at extremely 
low density, we are justified in neglecting the excluded volume contribution to the total 
free energy. 

Phase separation of an electrolyte or of a two component plasma is the result of 
an electrostatic instability arising from the strong positional correlations between the 
oppositely charged ions. This mechanism is very different from the one driving the phase 
separation in systems dominated by the short ranged isotropic forces. In that case the 
thermodynamic instability is a consequence of the competition between the interparticle 
attraction and the hardcore repulsion. 

The reduced temperature can be written as T* = a/XB, where Xb = (f/kBTe 
is the Bjerrum length. For water at room temperature ~ 7 A. This means that 
one would need ions of size less than 0.4 A, in order to observe phase separation at 
room temperature. This is clearly impossible since the minimum hydrated ionic size 
is about 2 — 4 A. Therefore, in order to see phase separation, one is required to look 
for solutions with As on the order of 40 A or more. For water such large values of 
Xb correspond to temperatures well below the freezing. An alternative is to work with 
organic solvents which have dielectric constants significantly lower than water. This 
was the strategy adopted by K.S. Pitzer in his studies of ionic criticality |5^, Q [5^ . 
Pitzer used liquid salt triethyl-n-hexylammonium triethyl-n-hexylboride (A^2226-B2226) 
in the diphenyl ether. With this he was able to observe the critical point at room 
temperature. Pitzer's work has provoked a lot of stimulating controversy because his 
measurements suggested that the Coulombic criticality belonged to a new universality 
class . At first sight this might not seem very surprising, after all the Coulomb force is 
extremely long ranged. On further reflection the situation is not so clear. Although the 
bare interaction potential between any two ions is long ranged, inside the electrolyte 
solution it is screened by the surrounding particles, as is seen from Eq. (0). The 
effective interaction potential, therefore, is short ranged, which should place the ionic 
criticality firmly in the Ising universality class. In fact all the theoretical arguments 
lead to this conclusion, which seems to be contradicted by the Pitzer's experiments. 
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In principle, it is possible that one has to be very close to the critical point before 
the Ising behavior becomes apparent. However, even this conclusion is hard to justify 
theoretically. Estimates of the Ginzburg criterion suggest that the width of the critical 
region for the Coulombic criticality should be comparable to that of systems with short 
ranged isotropic interactions The situation remains unclear. 

An alternative to working with electrolyte solutions is to study molten salts, which 
are classical two component plasmas. In this case the dielectric constant can be taken to 
be that of vacuum, and ions are no longer hydrated. The reduced critical temperature 
T* = 1/16 and the characteristic ionic diameter of about 2 A, imply that at criticality 
\b ^ 30, which means that the critical point for a molten salt is located at about SOOOi^'. 
It is, indeed, very hard to study critical phenomena at such high temperatures! It seems, 
therefore, that we are stuck with the low dielectric solvents. An alternative is the 
computer simulations, which are becoming sufficiently accurate to allow measurements 
of the critical exponents, at least for symmetric 1:1 electrolytes. Indeed the most recent 
simulations suggest that the Coulombic criticality belongs to the Ising universality 



class [57 



2.2. The Bjerrum association 

The DH theory presented in the previous section was based on the linearization of the 
Poisson-Boltzmann equation. In view of the strong screening and the rapid decrease of 
the electrostatic potential away from the central ion, such a linearization can be justified 
at intermediate and long distances. It is clear, however, that the linearization strongly 
diminishes the weight of configurations in which two oppositely charged ions are in a 
close proximity. Linearization underestimates the strength of electrostatic correlations 
which result in dipole-like structures. At low reduced temperatures characteristic of 
the critical point, these configurations should be quite important and must be taken 
into account. One way of doing this, while preserving the linearity of the theory, is to 
postulate the existence of dipoles with concentration governed by the law of mass action. 
In the leading-order approximation the dipoles can be treated as ideal non-interacting 



specie |5^, The total number of particles N = pV is then subdivided into 

monopoles A^^i = piV and dipoles N2 = P2V. The particle conservation requires that, 
= A''i + 2A^2- The free energy of the mixture is F = Ff"* + ^2^"* + F*^', where F^' and 
jpent ^YiQ entropic and the electrostatic free energies of the monopoles, given by the 
Eqs. (|1^) and ([T5|), but with N ^ Ni and p — ^ pi. The entropic free energy of dipoles 
is, 

PF^' = N2 ln[p2AVC2] - N2 , (24) 
where the internal partition function of a dipole is, 

C2{R) = 47r rr^dr exp ( ^] . (25) 



er ^ 

At low temperatures, the precise value of the cutoff R at which the two ions can be 
considered to be associated is not very important. Following the original suggestion 
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of Bjerrum|Q we can take this value to be the inflection point of the integral in 
Eq. (p5D, Rbj = Ab/2. This choice corresponds to the minimum of integrand in Eq. (|25|), 
which in turn can be interpreted as the probability of finding two oppositely charged 
ions at the separation r. The minimum then correspond to a liminal between bound 
and unbound configurations. A much more careful analysis of the dipolar partition 
function has been carried out by Falkenhagen and Ebeling based on the resummed 
virial expansion They found that that the low temperature expansion of the 

Bjerrum equilibrium constant is identical to the equilibrium constant which can be 
constructed on the basis of the resummed virial expansion. Since we are interested in 
the low temperature regime where the critical point is located, the Bjerrum equilibrium 
constant, (2 = (2(^3 j), will be sufficient. 

It is important to keep in mind that at this level of approximation the electrostatic 
free energy Ff"* is only a function of the density of free unassociated ions pi, since 
the dipoles are treated as ideal non-interacting specie. The concentration of dipoles is 
obtained from the law of mass action 

/i2 = yU+ + /i- , (26) 

where the chemical potential of a specie s is 
OF 

(27) 



V 



Substituting the expression for the total free energy into the law of mass action leads to 
P2 = \plC2e''''\ (28) 

where the excess chemical potential is /i*^^ = df^^/dpi. The critical point can be 
located from the study of the convexity of the total free energy as a function of ion 
concentration p. There is, however, a simpler way We observe that at Bjerrum 
level of approximation, dipoles are ideal non-interacting specie. This means that 
they are only present as spectators and do not interact with the monopoles in any 
way. This implies that only the monopoles can drive the phase separation. Thus, 
at the critical point the temperature must still be T* = 1/16 and the density of 
monopoles must still remain pl^ = I/GAtt, as in the case of the pure DH theory. 
The corresponding density of dipoles at criticality is then given by Eq. (pS]), with 
T* = 1/16 = 0.0625 and pl^ = I/GAtc = 0.00497. We find that at the critical point 
the density of dipoles is ~ 0.02. In the vicinity of the critical point there are many 
more dipoles than monopoles, P2c/ Pic ~ 4. Within the Bjerrum approximation the 
non-linear correlations, in the form of dipoles, do not affect the critical temperature, 
but strongly modify the critical density, p* = pl^ + 2p2c = 0.045. In spite of the 
crudeness of approximations, the location of the critical point agrees reasonably well 



with the Monte Carlo simulations [61, 62, 33|, T* = 0.051 and p* = 0.079. The 



coexistence curve, however, is found to have an unrealistic "banana" shape ||^. To 
correct this deficiency one must go beyond the "ideal" dipole approximation and allow 
for the dipole-ion interaction P, ||. Most of the fundamental physics of electrostatic 
correlations, however, is already captured at the level of the Bjerrum approximation. 
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3. Two-dimensional plasma and the Kosterlitz-Thouless transition 

An electrostatic system which over the years has attracted much attention is the two- 
dimensional plasma of positive and negative ions interacting by a logarithmic potential. 
The great interest in 2d plasma is due to the fact that various important physical systems 
can be mapped directly onto it. Examples include superfiuid ^He films, two-dimensional 
crystalline solids, and XY magnets ||6^. Although a continuous symmetry can not be 



broken in two dimensions [^, if the Hamiltonian of a system is invariant under an 
Abelian group, a finite temperature phase transition is possible. This transition occurs 
as the result of unbinding of the topological defects or "charges" . The defect-mediated 
phase transitions belong to the universality class of a two-dimensional plasma. 

Thirty years ago Kosterlitz and Thouless (KT) have presented a renormalization 
group study of the 2d plasma [^. They concluded that at sufficiently low temperature, 
the 2d plasma becomes an insulator. All the positive and negative ions pair-up forming 
dipoles. The metal-insulator transition was found to be of an infinite order, characterized 
by an essential singularities in the thermodynamic functions. The KT analysis, however, 
was restricted to the low ionic densities and it is not clear what happens when the 
concentration of charged particles is increased. It is tempting to apply to the 2d plasmas 



analysis similar to the one presented earlier for the 3d electrolytes . 

We shall, then, study a fiuid of disks with diameter a and charge ±g. The solvent 
is a uniform medium of dielectric constant e. The bare interaction potential for two ions 
(i, j) separated by distance r is 

= -Miin(r/a) . (29) 
e 

As in the case of 3d electrolyte, the mean-field contribution to the electrostatic free 
energy is zero, and all the important physics comes from the electrostatic correlations. 
We shall account separately for the long-ranged and the short-ranged correlations. The 
short ranged correlations lead to the formation of dipolar pairs of density p2, while the 
long ranged correlations produce the screening. As in the case of 3d electrolyte, the 
total density of hard discs is divided between the dipoles and the monopoles, so that 
p = Pi + 2p2- 

To calculate the electrostatic free energy we fix one ion and study the distribution 
of other particles around it. It is important to recall that the Poisson equation in 2d, 

V20(r) = -— p,(r) , (30) 



differs from the one in 3d by the normalization factor |67|, 27r has replaced the Air of 

the 3d Poisson equation. 

As before, we shall approximate the potential of mean force by the electrostatic 

potential and then linearize the Boltzmann factor. Linearization is compensated by the 

allowance for dipolar formation. 

The electrostatic potential for distances r < a satisfies the Laplace equation 
= 0, while for r > a the potential satisfies the Helmholtz equation = 
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with 

/2vrg2pi /on 

This equations can be easily integrated yielding the electrostatic potential, 

<P<{r) = — ln(r/a) + — — - for r <a, (32) 

e e KaKi[Ka) 

and 

0>(r)= for r>a, (33) 

e KaKi[Ka) 

where Ky{x) are the modified Bessel functions of order v. For large distances, the 
electrostatic potential decays exponentially fast. Just as in three dimensions, the 
electrostatic interactions are screened inside the two dimensional plasma. 



Eq. (^) consists of two terms, the potential produced by the fixed ion and the induced 
potential due to ionic atmosphere, 

* = . (35) 

e KaKi[Ka) 

Given the induced potential, the electrostatic free energy can be obtained using the 
familiar Debye charging process, Eq. (|12[). We find the electrostatic free energy density 
of a 2d plasma to be 

/3r' = ^\n[KaKi{Ka)] . (36) 

At the Bjerrum level of approximation, the electrostatic free energy depends only on 
the density of monopoles. The total free energy density is the / = /f"* + /|"* + /^', 
where 

/3/r* = Pi ln[piAV2] - Pi (37) 

and 

/3/r* = P2ln[p2A7C2]-p2. (38) 
The internal partition function for a 2d dipole is 

C2{R) = 27r / rdr exp In I — 

Ja L e Va. 

It is convenient to define the reduced temperature and density as T* = ksTe/q^ and 
p* = pa?. We note that for low temperatures, T* < 1/2, the integral in Eq (^) 
converges uniformly as R —>■ oo. In this regime it is possible, therefore, to define the 
internal partition function of dipole as 

^ , , 27ia^T* 

(2 = C2(oo) = . (40) 



(39) 
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Figure 3. Phase diagram for tlie two dimensional plasma within the Debye-Hiickel- 
Bjerrum approximation. We expect the fluctuations to renormalize the Kosterlitz- 
Thouless line, shifting it from its horizontal position and making it density dependent. 
The topology of the phase diagram should, however, remain the same. 



The thermodynamic equilibrium requires that for fixed volume and number of particles 
the Helmholtz free energy be minimum. This is equivalent to the law of mass action 
Eq. (p6D, which upon the substitution of free energy simplifies to Eq. (|28|). In the limit 
of small concentrations, the excess chemical potential can be expanded in powers of pi 
yielding 

^3^^'' = -^bE + H'^a/2)], (41) 

where •je is the Euler constant. Substituting Eq. (^T]) into Eq. p8|), we find that the 
concentration of dipoles in the limit p ^ scales as 

P2 ~ P?^*^ , (42) 

where 

(^in = 2 - ^ . (43) 

For T* < 1/4, the exponent 9{T*) < 0, and in the limit pi ^ the law of mass action 
can not be satisfied. This means that in the temperature density plane (T*,p*), for 
sufficiently small densities, the line T* = 1/4 corresponds to the critical locus of metal- 
insulator transitions. Below this line, and for sufficiently small ionic concentrations, no 
free monopoles can exist. All the ions are paired up into neutral dipolar pairs. The 
critical line terminates at the tricritical point located at T^j, = 1/4 and 

Ptri = -7, — - 0.003954 . (44) 

For T* < 1/4 and p* > pj'^j there is a phase separation between an insulating vapor 
and a conducting liquid phases. Fig ^ As the critical line is approached from high 
temperatures, the Debye length diverges as 

= ~ e'=('')/*^ , (45) 
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where 



cip) = ]\n(P^] , (46) 



and 



4 V P 

< = ^^. (47) 

The critical exponent is z/ = 1. The Kosterlitz-Thouless (KT) renormalization group 
theory |^ predict the same behavior for C,d except that u = 1/2. The KT theory, 



however, leaves unanswered the question of what happens to the metal-insulator 
transition for higher plasma concentrations. The current theory, on the other hand, 
shows that the critical line terminates in a tricritical point, after which the metal- 
insulator transition becomes first order |]6^. This topology is also consistent with the 



findings of Monte Carlo simulations [^, ^ . A more sophisticated theory introduced by 
Minhagen ||70|, leads to a very similar phase diagram, except that the tricritical point 
is replaced by a critical end-point. 

We see that the electrostatic correlations are even more important in 2d than in 
3d. While in three dimensions the electrolyte phase separates into the coexisting liquid 
and gas phases, both of which contain monopoles and dipoles, in two dimensions the 
low density vapor phase does not contain any free charges and is an insulator. 

4. The one component plasma 

The one component plasma {OCP) is probably the simplest model of a Coulomb system. 
It consists of point ions, all of the same sign, inside a rigid neutralizing background of 
opposite charge. In spite of its simplicity the model is relevant for many physical systems. 
Some examples are the interiors of stars, liquid metals, and magnetically confined 
electrons. As an approximation, OCP is particularly important since it provides a 
framework in which ionic correlations can be calculated. Thus, the electrostatic free 
energy of a homogeneous OCP can be used to account for the counterion correlations 
in colloidal suspensions and polyelectrolyte solutions. The OCP has been extensively 
studied over the years. A review of the subject, which still remains very actual, has 



been presented by Bans and Hansen some twenty years ago [fTlH . 

In the spirit of the current work we shall, however, confine our attention to simple 
analytical theories of the OCP ||72|, [73[. Our model consists of ions each carrying 



charge q inside a uniform neutralizing background of dielectric constant e and volume 
V. The average ionic density is p = N/V. The mean electrostatic potential inside the 
OCP is zero and the free energy is, once again, entirely due to positional correlations 
between the ions. The advantage of working with the OCP is that it contains only two 
independent length scales: the average separation between the particles d = (AiTp/3)~^^^ 
and the Bjerrum length Xb = q^/kBTe. There is only one dimensionless parameter on 
which all the thermodynamic quantities depend, F = Xs/d. This is quite distinct from 
the electrolytes and two component plasmas, for which besides d and there is a third 
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length scale — the ionic diameter a — and the thermodynamics is parametrized by two 
dimensionless quantities T* and p*. For electrolytes the limit a —>■ does not exist, since 
for point particles the free energy can be lowered indefinitely by collapsing the system 
into point-like dipolar pairs. 

The calculation of free energy of the OCP proceeds along the lines of the Debye- 
Hiickel theory. We fix one ion and study the distribution of other ions around it. The 
electrostatic potential satisfies the Poisson equation, Eq. (|), with the charge density 
given by 

Pq{r) = qpg{r) - qp . (48) 

The correlation function can be expressed in terms of the potential of mean force w{r) 
as 

g{r) = e-^""^'^ . (49) 

In the spirit of the Debye-Hiickel theory we replace w{r) by q(j){r). This approximation 
entails neglect of electrostatic correlations inside the ionic cloud which surrounds the 
central ion. The approximation should be quite good as long as F is not too large. 
As the next step, linearization of the Boltzmann factor leads to the charge density of 
particularly simple form, 

P,(r) = > (50) 

where = inPq'^p/e. Substituting this into the Poisson equation we, once again, 
find the familiar Helmholtz equation (^. This can be easily integrated yielding the 
electrostatic potential of Yukawa form, 

0(r) = ^ . (51) 
er 

The charge density must be bounded from below pg > —qp. However, considering 
Eqs. ( |50[) and (|5lD , it is evident that this condition is violated for sufficiently small 
separations from the central ion. Something must have gone seriously wrong. It is 
easy to trace the problem to the linearization of the Boltzmann factor. Clearly at short 
distances linearization is not justified since the strong electrostatic repulsion between the 
ions results in a very large electrostatic energy. It appears that in order to understand 
the physics of the OCP one has to solve the full non-linear Poisson-Boltzmann equation, 
a task which cannot be performed analytically. Fortunately not everything is lost. A 
way out was suggested by Nordholm, who noted that the strong repulsion between like- 
charged ions results in an effective hole surrounding the central ion []72|. Very few ions 
can penetrate inside the hole since this costs them too much electrostatic energy. Inside 
this correlation hole, r < h, the electrostatic potential satisfies the Poisson equation 
(I) with a uniform background charge density, Pg ~ —qp- In the outside region, where 
the electrostatic interaction is much weaker, the linearization of the Boltzmann factor 
is justified and the potential satisfies the Helmholtz equation. Therefore, for r > h the 
electrostatic potential still has a Yukawa form, but with an undetermined prefactor, 

P,ir) = ^— . (52) 
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The coefficient A can be obtained from the condition of continuity of Pq{r) across r = h. 
We find 

A = -qphe^^ . (53) 

The charge neutrality requires that 
47r f°° 

q qph^ + A-K I r'^drpg{r) = . (54) 

3 Jh 

The integral can be performed explicitly and the resulting equation solved to determine 
the size of the correlation hole 



h = d[uj{T) - 1]/V3r , (55) 

where 

u{r) = [1 + {srf/Y' ■ (56) 

The size of the correlation hole is a monotonically increasing function of the coupling 
strength. For high temperatures, small couplings, ~ A^. This is exactly what one 
might have expected, since on this length scale the electrostatic repulsion becomes 
comparable to the thermal energy. When the temperature is lowered the kinetic energy 
diminishes and the particles are strongly scattered by the electrostatic repulsion, causing 
increase in the size of the correlation hole. For very low temperatures, large F, the size 
of a correlation hole is equal to the average spacing between the particles, h ^ d. This 
is, once again, the correct limiting behavior since at low temperatures the ions tend to 
keep as far away as possible from their neighbors. For r > h the electrostatic potential 
follows directly from Eq. (0), 

^ , , Anqphe-''^''-''^ 



(57) 



For r < h the Poisson equation with a uniform background charge Pg = —qp must be 
solved. The solution is easily found to be 

^,(,.) = ^ + H^ + ^. (58) 

er 3e 

The induced potential ip can be determined from the condition of continuity of 
electrostatic potential 4>>{h) = 4><{h), which reduces to 

^ = -^{[l + (3r)3/T/^-l}. (59) 
2q 

The Debye charging process Eq. (|12D, yields the electrostatic free energy per particle 



N 

1 
4 



2 2tt /uj^ + iu + l\ 2 _i/2cu + l 
cj^ H ^ + In ^ tan ^ 



(60) 



3v^ V 3 y V 
where u!{r) is given by Eq. (|5B|). In the limit of high temperatures, F — 0, Eq. 
reduces to the Debye limiting law, Eq. (|14|). Furthermore, the free energy agrees with 
the Monte Carlo simulations with an error of less than 10% over a wide range of coupling 
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strengths, < F < 80. This suggests that the inclusion of a correlation hole into 
the Debye-Hiickel theory captures most of the essential physics of the one component 
plasma. 

Some caution must be taken when using the OCP to model real physical systems. 



For F > 3, the isothermal compressibility and pressure of the OCP become negative ||71|| . 
This is a consequence of treating the background as a rigid entity and neglecting its 
pressure. How this can be corrected, in practice, depends on the kind of problem that 
one wants to study. If one wants to use the OCP to model dense ionized matter, the 
suitable background is the degenerate electron gas. When the free energy of background 
is added to the OCP the pressure and the compressibility remain non- negative for all 
values of F. An alternative approach was suggested by Weeks, who defined the OCP 
as the classical "dense-point limit" of a two component plasma [Q. In this limit, 
the background is treated as an infinitely dense cloud of point particles each carrying 
an infinitesimal charge e, so that epback remains constant. Presence of such background 
does not affect the electrostatics of the OCP, but regularizes its pressure and isothermal 
compressibility, making them non-negative. An interesting byproduct of this analysis 
is the conclusion that the freezing of the OCP occurs without any change in density, 
i.e. the volume per particle in the fluid and the solid phases is the same WM. The 



freezing transition happens at F ~ 180 and the resulting solid phase has the BCC 



structure |7^, 171], |76| 



4.1. Confined one component plasma 

In 1971 Crandall and Williams suggested that electrons trapped on the surface of liquid 
helium ^He can crystallize, forming a two dimensional Wigner crystal [^. Eight 



years later this order- disorder transition was observed experimentally by Grimes and 



Adams ||7^. In this system electrons obey the classical mechanics, since the Fermi 
energy is much smaller than fc^T. Similar crystallization can occur in the inversion 
layer near the surface of a semiconductor, however, in this case the quantum effects are 
important and the electrons form a degenerate quantum gas [[79| . 

The trapped electrons above the liquid ^He can be modeled as a confined quasi- 
two-dimensional plasma of particles interacting by 1/r potential. This model is also 
appropriate for the study of correlations between the condensed counterions on the 
surface of colloidal particles. 

The average spacing between the confined electrons is = (vra)^^/^, where a is 
the average surface density, a = N/A. The dimensionless quantity parameterizing the 
strength of electrostatic interactions is F = q^/eksTd. For an infinitesimally thin layer 
separating two mediums of dielectric constants ei and 62, the important parameter is 
the average dielectric constant e = (ei -|- e2)/2. It has been observed in computer 



simulations [£0| that the 2d OCP crystallizes into triangular Wigner crystal for F > 130. 
This value is also in a close agreement with the experiments of Grimes and Adams. 
We can gain much insight into thermodynamics of 2d OCP using the, now familiar. 



20 



Electrostatic correlations: from Plasma to Biology 



Debye-Hiickel theory. Our model consists of a plasma of point particles of charge q and 
of a neutralizing background, confined to an interface located at z = between the two 
dielectric half-spaces. For 2; < the dielectric constant is ei and for z > the dielectric 
constant is €2. Since the half-spaces do not contain any free charges, the electrostatic 
potential everywhere satisfies the Laplace equation = 0. The electrostatic free 
energy is obtained by fixing one particle and calculating the induced potential resulting 
from the redistributions of other ions in the z = plane. It is convenient to adopt the 
cylindrical coordinate system, {g, if, z), so that the fixed ion is located at ^ = 0, z = 0. 
Using the azimuthal symmetry and the fact that the electrostatic potential vanishes at 
infinity, the solution to Laplace equation can be written as fSlI] 



0i(^,z)= / Ai{k)Jo{kg)e''' dk for z<0, (61) 

and 







POO 

02(^,^)= / A2{k)jQ{kQ)e-^' dk for z>0 , (62) 
Jo 

where Jq{x) is the Bessel function of order zero. 

The functions Ai{k) and A2{k) can be determined from the boundary conditions, 
which are: continuity of the electrostatic potential, 

02(^?,o) = 0l(^?,o) , (63) 

and discontinuity of the displacement field across the 2; = plane. The discontinuity 
results from the inhomogeneous distribution of interfacial charge induced by the fixed 
ion, 

[e2Ei2{Q,z)-eiEi{Q,z)]-n = A^aq{Q) . (64) 

From charge neutrality the average interfacial charge is zero so that crq(^) is the result 
of ionic correlations, 

= ^ -qa + gae''^^'^^^'") . (65) 
27ig 



The first term of Eq. (|65| ) is the surface charge density of the fixed ion, the second term is 
due to the uniform negative background, while the last term is the surface charge density 
of ions confined to the interface. We have, once again, approximated the potential of 
mean force by the electrostatic potential. In the spirit of Debye-Hiickel theory we shall 
now linearize the Boltzmann factor. The surface charge density becomes 

q6{g) e0(g,O) 

= -TT— - IT-T , (66) 

zng zhAgc 



where 



Agc = 7^ — 5- (67) 



27rg2cr 

is the Gouy-Chapman length. 
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The continuity of electrostatic potential requires that Ai{k) = A2{k). Substituting 
Eqs.([61|) and (p2| ) into Eq.(p^) and using Eq.(p6D, we find the electrostatic potential 
over the full range — oo < 2; < oo to be 

y / T fu.\ ^-k\z\ 



z) = ^ Mkg) e-'^^^dk . (68) 

For z = the integral can be performed explicitly yielding 

<p{q,0) = , (69) 



where the functions Tj^{x) are defined as |p2 

Tuix) = 1 - ^^[H,{x) - N,{x)] , (70) 

with Hi,{x) and Ny{x) being the Struve and the Bessel functions of order z/, respectively. 
For large values of x, tq ~ 1/a;^, so that asymptotically, 

<^(^,0)-%. (71) 

We conclude that in the case of a confined plasma there is no exponential screening, 
instead the electrostatic potential is purely algebraic and has the form of a dipole-dipole 
interaction |7T . 



There is a well known argument in the condensed matter physics going back all 
the way to Bloch Peierls and Landau in the 1930's, which states that a 
continuous symmetry can not be broken in two dimensions. This means that there can 
not exist a true two dimensional crystalline order, since it requires breaking translational 
symmetry. The argument was made rigorous by Mermin, who proved it for particles 
interacting by short-ranged potentials [^]. It is quite simple to see how this conclusion 
arises. Suppose that there is a 2d crystal, one can then calculate the mean-square 
displacement 5^ of one particle from its equilibrium position due to thermal fiuctuations. 
It is found that 5^ ~ TlnL, where L is the characteristic crystal size. For L 00, 
the mean square displacement diverges for any finite temperature, implying that in 
thermodynamic limit a 2d crystal is unstable to thermal fiuctuations. Although there 
is no true long-range order in two dimensions for systems with short-range forces, there 
exists a pseudo-long-range order characterized by an algebraically decaying correlation 
functions. It is not clear, however, to what extent this conclusion applies to the 2d 
OCP, whose particles interact by a long-ranged 1/r potential. Certainly in this case 
Mermin's proof is no longer valid. However, since the effective interaction potential 
inside a 2d OCP decays as 1/r^, which is short-ranged in two dimensions, suggests that 
there should not be any long-range order. Whether there is a true long-range order or 
a pseudo- long-range order for a 2d OCP remains uncertain. Simulations find that for 
r ^ 130 there is a crystallization transition. It is, however, difficult to say whether 
the crystalline state has a true long-range order or a pseudo-long-range order [|0|. It IS 
also unclear if the transition is of first order or continuous, belonging to the Kosterlitz- 
Thouless universality class Existence of the thermodynamic limit for confined 
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2d plasmas can also be attributed to the effective renormalization of the interaction 
potential from a non-integrable 1/r to integrable (in two dimensions) form. 

The Helmholtz free energy of a 2d plasma can be obtained directly from Eq. (|69|). 
We need to know the induced potential felt by the central ion due to other particles. In 
the limit ^? — *• 0, the electrostatic potential reduces to 

0(^,0) ^^ + -^ln(f5/2AGC7)- (72) 
eXcc 

The first term of this expression is the potential produced by the central ion, while the 
second term is the induced potential felt by the fixed ion. We note that the induced 
potential is actually divergent in the limit ^ ^ 0. This is the consequence of the failure 
of linearization of the Poisson-Boltzmann equation. This deficiency can be corrected in 
the same way as was done for the 3d OCP, by introducing a correlation hole of radius h 
around each particle. Unfortunately in the present geometry this leads to calculations 
which are no longer tractable analytically. From our study of the 3d OCP we can, 
however, make some reasonable approximations. In the limit of high temperatures, 
small r, the size of the correlation hole should be such that the electrostatic and the 
thermal energies become approximately equal, 

f^^ksT. (73) 

This means that h X^. We can use this value as the short-distance cutoff in the 
calculation of free energy. The induced potential then becomes 

^~^ln(AB/2AGc). (74) 

The free energy is obtained through the usual Debye charging process, Eq. (|12|). 
Recalling that AB(Ag) = A^As(g) and \Gc{.^<i) = Agc(q')/A^, where A is the charging 
parameter, in the limit of high temperatures F — > 0, the reduced free energy per particle 
is found to be 

^^F^MF). (75) 

Eq. (^) is precisely the leading order term of the resumed virial expansion obtained by 
Totsuji 

For low temperatures, the OCP crystallizes into a triangular lattice. The Madelung 
energy of this lattice is, 

^ = -1.106103F . (76) 

This equation provides a surprisingly good fit not only for the free energy of solid, but 
also for the free energy of fluid at sufficiently high values of F. Comparing to the results 
of the Monte Carlo simulations we find that for F = 5 the error accrued from using 
Eq. ( [7B| ) to calculate the total electrostatic free energy is about 30%. For F = 20 this 
error drops to 11% and for F = 50 it goes down to 6%. Recalling that the crystallization 
transition occurs at F ^ 130, we see that the Eq. (|7^) works well into the fluid phase. It 
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is reasonable, therefore, to approximate the electrostatic free energy of a fluid for F > 5 
by 

-j^ = -1.106103r . (77) 

The reason why the electrostatic free energy of a fluid is so well approximated by the 
free energy of the crystal, is a consequence of strong electrostatic correlations. 

5. Asymmetric systems 

Up to now we have considered only symmetric plasmas and electrolytes. In practice, 
however, it is unlikely that both cations and anions have exactly the same size and 
magnitude of charge. It is, therefore, important to explore the thermodynamics of a 
general Z : 1 electrolyte in which cations have charge Zq and diameter Oc, while anions 
have charge —q and diameter a^. Unfortunately, as soon as the asymmetry is introduced, 
the internal inconsistency enters into the Poisson-Boltzmann equation ^0 . Recall that 



the cation-anion correlation function can be expressed in terms of the potential of mean 
force w+_ 

g^_(r) = e-^^^+^W . (78) 

The w+_(r) is the work needed to bring cation and anion from infinity to separation r 
inside an electrolyte. Clearly this work is invariant under the permutation of particle 
labels (r) = w_+(r). This means that 

9+-{r) = 9-+{r) . (79) 

The Poisson-Boltzmann equation, which serves as the basis, for the Debye-Huckel theory, 
approximates the potential of mean force by w+_(r) = g_0+(r). The self consistency 
condition, Eq. (^), then requires that 

g+0_(r) = g„0+(r) . (80) 

Because of the non-linear nature of the PB equation this condition can not be satisfied 
except for symmetric electrolytes. The linearization prescription intrinsic to the Debye- 
Hiickel theory allows Eq. (|80|) to hold for ions of different valence, but with the same 
ionic diameter, = a^. 

We see that as soon as the symmetry between the cations and anions is broken the 
physics and the mathematics of the problem becomes significantly more complex. In 
the limit of very large asymmetries, Z ^ oo and Oc ^ a new simplification, however, 
enters into the game. 

5.1. Colloidal suspensions 

A typical colloidal suspension often studied experimentally consists of polystyrene 
sulphonate spheres of diameter lOnm — l/im and 10^ — 10^ ionizable surface groups. 
Because of the large surface charge, the colloidal particles tend to repel each other. 
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forming crystals, even at fairly low volume fraction of less than 10%. Using the periodic 
structure of the lattice, the thermodynamics of a colloidal crystal can be studied fairly 
straightforwardly. Each colloidal particle can be thought to be confined to a Wigner- 
Seitz (WS) polyhedral cell. A further approximation replaces the polyhedral WS cell 
by a sphere [p2| . 



5.2. Colloidal lattices 

We shall model the colloidal particles as hard spheres of radius a carrying Z ionizable 
groups of charge —q distributed uniformly on the surface. The counterions will be 
idealized as point particles of charge +q. The suspension of Np = ppV polyions and 
Nc = ZNp = pcV counterions is confined to a volume V. As usual, the solvent will be 
treated as a uniform continuum of dielectric constant e. For sufficiently large polyion 
concentrations colloidal suspension crystallizes. Using the lattice symmetry, we restrict 
our attention to one colloidal particle and its counterions inside a spherical WS cell of 
radius R such that 

Pp = 1^ • (81) 

Using the statistical mechanics it is possible to show that the osmotic pressure inside a 
cell is proportional to the concentration of counterions at the cell boundary |]88| , 



(3P = p,{R) . (82) 

The thermodynamics of a crystalline colloidal suspension now reduces to the calculation 
of the distribution of counterions inside a WS cell. This can be done using a simple 
mean-field picture. The electrostatic potential inside a WS cell satisfies the Poisson 
equation (Q) with the counterion charge density approximated by the normalized 
spherically symmetric Boltzmann distribution, 

e-/3q</'(0 

p,(r) = ZqNp — - . (83) 

J a 

The non-linear Poisson-Boltzmann equation can be solved numerically to yield the 
electrostatic potential and the distribution of counterions inside the cell. In practice it 
is more convenient to work with the electric field 

E(r) = -V0(r) . (84) 

The Poisson equation can then be rewritten as 

V-E(r) = ^[p,(r)+g,(r)], (85) 

where ^^(r) is the polyion charge density, 

?p(r) = -£^5(|r|-a)]. (86) 
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To simplify the calculations we have uniformly smeared the charge of the polyion over 
its surface. Integrating both sides of Eq. (p5| ) and taking advantage of the divergence 
theorem, the electric field at distance r from the polyion is 

E{r) = --^ [Zq - air)] , (87) 



J |r'|<r 



where 

' \r'\<r 

is the counterion charge inside a sphere of radius r centered on the colloidal particle. 
Using the gauge in which 0(a) = the electrostatic potential is 

0(r) = - / drE{r) (89) 

J a 

and the Poisson-Boltzmann equation reduces to an integral equation for the electric 
field. Note that Eq. (|87| ) naturally incorporates the boundary conditions 



E(a) = (90) 

and 

E(R) = . (91) 

Eq. ( p7D can be solved iteratively to yield the counterion density profile from which all 
other thermodynamic functions are straightforwardly determined. 

For aqueous colloidal lattices with monovalent counterions, the Poisson-Boltzmann 
equation is in excellent agreement with the experiments and simulation. The PB 
equation, however, does not account for the correlations between the counterions and 
breaks down for low dielectric solvents or for aqueous suspensions with multivalent ions. 
Fortunately, in the case of colloidal lattices, it is fairly straightforward to account for 
these effects using the density functional theory. 

5.3. Density functional theory {DFT) for colloidal lattices 

We shall now construct the Helmholtz free energy functional for a Wigner-Seitz cell 
of a colloidal lattice. The Helmholtz free energy is a functional of the average local 
counterion concentration 

Pe(r) = (5^5(r-r.)), (92) 

1=1 

where the brackets denote a Boltzmann average over all the particle positions. The free 
energy consists of electrostatic and entropic contributions. The entropic contribution is 
simply that of an inhomogeneous ideal gas, 

/3F-*[p,(r)] = I rf3rp,(r){ln[p,(r)A3] - 1} . (93) 
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The electrostatic contribution is the result of Coulomb interactions between the 
counterions and the polyion, as well as the self-energy of the ionic cloud, 

i^iPc(r)] = g / dhdh' [ dhdh' . (94) 

J e\r — r'l 2 J e|r — r'| 

Eq. ( |94D is the mean-field approximation for the total electrostatic free energy. It does 
not account for the electrostatic correlations between the counterions surrounding the 
colloidal particle. Clearly, if there is a counterion present at position r, there is a reduced 
probability of finding another counterion in its vicinity. This information is not included 
in Eq. (P^. One can attempt to account for the electrostatic correlations using the local 
density approximation (LDA), 

/3F-'^[p,(r)] = I dhp,{r)U.[Pc{r)] . (95) 

The correlational free energy can be approximated by the free energy of a three 
dimensional one component plasma Eq. (|60|) , 

fcor[pc{r)] ~ focp[pc{r)] = 

1 

4 



io^ H T= + In ^ tan ^ ' 



(96) 



with 



3^3 V 3 y V3 \ Vs 

u;[p^{r)] = [1 + {3T[p^{r)]fY' ' (97) 
where the local coupling strength is, 

r[p,(r)] = AB[47rp,(r)/3]V3. (98) 

Although Eq. (|96D was derived for a plasma in a uniform neutralizing background, it 
can also be used to approximate the correlational electrostatic free energy of counterions 
confined in a WS cell. The role of a neutralizing background being played by the 
confining electric field produced by the colloidal particle. 

The equilibrium charge distribution is determined from the minimization of the 
total Helmholtz free energy 

F = F^"* + F^' + F"'"' (99) 

subject to constraint of particle conservation 

jdhp,{r) = N,. (100) 

This is equivalent to minimization of the grand potential 

n = F - /ieiV, , (101) 

where the chemical potential of counterions pc is the Lagrange multiplier. Performing 
the calculation, leads to the equilibrium counterion density profile 

-/3g0(r)-/3/.'=-{r) 

Pc{r) = N, H , (102) 

Ja 
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where the excess chemical potential is, 

/i'^^(r) = . (103) 

We see that in the absence of correlations the density profile reduces to the Boltzmann 
distribution Eg . (|83|) . The counterion density, which enters into the expression for the 
excess chemical potential, can be expressed in terms of the electric field, 

Pc(r) = ^V-E, (104) 
which due to spherical symmetry simplifies to 

= . (105) 

47rgr^ or 



Substituting Eg . (11021) into Poisson eguation (|85D and using Egs. (|89|) and (|105D leads to 



an integro-differential eguation for the electric field, Eg.( 

Unfortunately, for large colloidal charges this eguation has no stable solutions [^, 
pO| , There is an unbounded increase in the concentration of counterions in 

the vicinity of the colloidal surface resulting from the failure of the local density 
approximation |Q. To overcome this difficulty Groot and |^ has suggested use 



of a weighted density approximation {WD A), which has proven guite successful in its 
application to other problems of the condensed matter physics [^, Within this 
approach the correlational free energy is given by 

F-lp,(r)] = I dhp,{r)Ur[Pc{r)] , (106) 

where Pc(i') is the average local density. 



p,(r) = j dh'w{\r-r'\)p,{r) . (107) 

The weight function w{r) can be determined from the thermodynamic reguirement that 

S^BF 6(r-r') 
, = -^YT^ - C^ilr - r'l) , (108) 

dpc{r)6pc{r'} pc{r) 

where C2(|r — r'|) is the direct correlation function. In particular this eguation must 
hold in the limit of a homogeneous OCP, the direct correlation function for which can 



be obtained using the theory developed in Section | |9l 

h 

and 



C2(r) = if r<h, (109) 



C^ir) = -— if r>h. (110) 
r 

Performing the calculation, we find that the weight function is well approximated by 0] 

where G(x) is the Heaviside step function, Q{x) = 1 for x > and 6(x) = for 
X < 0. The WD A is significantly more computationally demanding than the LDA. Its 
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advantage, however, is the numerical stability for all values of the colloidal charge. The 
results based on the WD A are in excellent agreement with the Monte Carlo simulations. 
The WD A gives us a good handle on the thermodynamics of colloidal lattices. At lower 
concentrations, when the crystalline structure has melted, the situation unfortunately 
is no longer so clear cut. In this case a simple picture based on the Wigner-Seitz cell is 
not sufficient and new methods must be developed |^, ^ . Unfortunately the standard 
techniques of the liquid state theory based on integral equations are powerless in the 
case of highly asymmetric colloidal systems. The field theoretic methods also fail when 
applied to this difficult problem. Fig. |l|. Furthermore, even the experimental situation is 
far from clear. Ise and coworkers claim to have seen stable clusters of colloidal particles 
in highly deionized colloidal suspensions. Tata et al. even report an observation of a full 
equilibrium vapor-liquid-like phase separation [Q. These experiments, however, have 
been challenged by Palberg and Wiirth, who demonstrated that the phase separation 
observed by Tata et al. was the result of non-equilibrium salt gradients produced by 
the ion exchange resin |9^, In the colloidal science community the possibility of a 
liquid-vapor phase separation in highly deionized colloidal suspensions has met with a 
large amount of scepticism. The usual argument against the phase transition is based 
on the Derjaguin-Landau-Verwey-Overbeek (DLVO) colloidal pair potential [^, |10CI| . 



It is easy to understand the nature of the DLVO potential based on the Debye- 
Hiickel theory. If the size of colloidal particles is shrunk to zero, a — > 0, then due to 
screening by counterions, the interaction energy between two "point" colloids would be 
of a Yukawa form, 

p-Kr 

Voir) = {Zqf^— , (112) 
er 

where the inverse Debye length is 



Now, consider the electrostatic potential outside the fixed colloidal particle of radius a 
and charge —Zq^ Eq. ( pUj ) 

0> = -^-^ , = TT— ^ • 114 

er (1 + x) 

Evidently the factor 9{na) accounts for the fact that screening starts only outside the 
cavity, r > a. We also can think of Eq. (|1 141) as the potential of a point particle with 
an effective charge Qp = Zq6{K,a). An advantage of this alternative point of view is 
that the interaction energy for two "point" particles is simply given by Eq. ( |112| ) with 
Zq Qp. This leads directly to the famous DLVO potential 

-ftr 

VDLVoir) = {ZqfQ\Kay^ . (115) 

er 

This potential is purely repulsive ||1(J1[ |, which naively suggests that a charged colloidal 
suspension is stable against a liquid-gas phase separation. Sogami and Ise, therefore, 
have argued that the DLVO potential must be incorrect, since it cannot account for the 
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inhomogeneities observed experimentally [ p.02|| . In its stead, they proposed a different 
interaction potential derived on the basis of the Gibbs free energy. The potential found 



by Sogami and Ise contains a minimum |102|, which implies that at short enough 



separations the two like-charged colloidal particles attract! What is most surprising is 
that the attraction appears even for monovalent counterions, i.e in the absence of strong 
correlations between the colloidal double layers. Furthermore, water is an incompressible 
fluid so that it is difficult to see how a change of paradigm from Helmholtz to the 
Gibbs free energy can lead to such a profound modification of the interaction potential. 
Inconsistency in the results based on the Helmholtz and the Gibbs free energies has 
been carefully reexamined by Overbeek, who has traced the discrepancy to a flaw in the 
Sogami and Ise's calculations [p.03| . 



It is important to stress that the repulsive two-body interactions do not, in general, 
preclude the possibility of a liquid-gas phase separation in a multicomponent fluid. 
In fact van Roij and Hansen found, within the linearized density functional theory, 
that it is possible for a colloidal suspension with polyions interacting by the repulsive 
DLVO potential to phase separate into coexisting liquid and gas phases ||2^. Before 



entering into the discussion of colloidal fluids it is, however, important to introduce a 
new fundamental concept — the colloidal charge renormalization. 

5.4- Charge Renormalization 

Although the non-linear Poisson-Boltzmann equation can not be solved analytically for 
a spherical geometry, the numerical solution indicates that the electrostatic potential 
far from colloidal particle saturates function of the bare colloidal charge ||22|| . 



This suggests that the thermodynamics of a highly charged colloidal systems can be 
based on a linearized PB equation but with the bare colloidal charge replaced by an 
effective renormalized charge. The original concept of colloidal charge renormalization 
is due to Alexander et al., but is well predated in the polyelectrolyte literature, where 
the phenomenon is known as the Manning counterion condensation ||104| , |105| , |106| . 



To understand better colloidal charge renormalization, let us first consider a 
uniformly charged plane at fixed potential ips inside a salt solution of concentration 
c. The electrostatic potential at distance x from the plane satisfies the PB equation, 

— = smh(/?g0) . (116) 

Since at the moment we are considering aqueous suspensions containing only monovalent 
ions, the electrostatic correlations are insignificant and the mean-field Poisson- 
Boltzmann approximation is sufficient. Multiplying both sides of Eq. (|116|) by d(j)/dx 
allows us to perform the first integration. Since the potential vanishes in the limit 
a; — > oo, we find 

lW{x)r = ^[cosh{Pq<P)-l] . (117) 
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The second integration yields ||107|| 



2kBT 1 + e-nanh(/3gV;3/4) 

(pix) = m- — — — -^-r , (118) 

^ q 1 - e-^^ tanh(/5gV^,/4) ' ^ ^ 

where the inverse Debye length is, 

K = ^/SttcXb ■ (119) 

In the limit of large surface potentials this expression simplifies to 

2kBT 1 + e--" 

'^(^) = ^l^ l_e-.x • (120) 
For separations from the plane larger than the Debye length, Eq. (|120|) becomes 

(f){x) = ^hlLe-- . (121) 
Q 

An important observation is that for large surface potentials, f3qips/4: ^ 1, the 
electrostatics away from the plane is completely insensitive to the surface charge density. 

Now, let us consider a highly charged colloidal particle of valence Z and radius 
a inside a symmetric 1 : 1 electrolyte of concentration c. The electrostatic potential 
at distance r from the center of a colloidal particle satisfies the PB equation (^). For 
distances r > a + C,d the electrostatic potential is small and the PB equation can be 
safely linearized leading to the Helmholtz equation (|^). This can be easily integrated 
yielding the electrostatic potential, 

-Kr 

0(r) = A , (122) 

r 

where A is the integration constant. To find its value, lets restrict our attention to 
suspensions in which the <^ a. In practice this is not a very strong restriction. 
For salt solutions at physiological concentrations ~ 8 A while the characteristic 
colloidal size is on the order of 1000 A. Even for solutions with very low salt content, 
in the mM range, the Debye length is on the order of 100 A. Under these conditions 
all the curvature effects associated with the spherical geometry of colloidal particle are 
effectively screened at separations a + ^D < ^ < 2a, and the electrostatic potential is well 



approximated by that of a uniformly charged plane, Eq. ( 121 ). Comparing Eqs. ( |121| ) 
and Eq. (|122|) the value of the integration constant follows directly, and the electrostatic 
potential at distance r > a + C,d from the center of colloidal particle is 

0(r) = — ^ . 123 

qr 

This is the asymptotic solution of the full non-linear PB equation for ^ 1. 
Comparing this to the solution of linearized Pi?, Eq. (114), it is evident that the two 



are identical as long as the hare colloidal charge is replaced by the renormalized charge. 
For highly charged particles, Eq. ( |123[ ) shows that the renormalized charge saturates 
at 

^;e1-^^^. (124) 
Ar 
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While the previous analysis was carried out for one colloidal particle inside an 
electrolyte solution, the concept of charge renormalization is quite general and can be 



applied to colloidal suspensions under various conditions ||109| , p.09| , |110| , |111|| . The 
difficulty of defining the effective charge for suspensions at non-zero concentrations 
resides in the complexity of accounting for the consequences of colloidal interactions. 
The standard practice is to to study one colloidal particle inside a spherical Wigner- 
Seitz cell whose radius is determined by the volume fraction of colloids ||2^. While this 
procedure is fully justified for colloidal lattices, its foundation is less certain for fiuidized 
suspensions. To find the renormalized charge one numerically solves the full non-linear 
PB equation and matches the electrostatic potential to the solution of the linearized 
equation at the cell boundary. Alternatively, the osmotic pressures inside the WS cell 
calculated using the non-linear and linear equations are matched in order to define the 
effective charge. One should remember, however, that while at the level of non-linear 
PB equation the osmotic pressure is directly proportional to the concentration of ions 
at the cell boundary, Eq. (^2]), this is not the case for the linearized PB equation. The 
various procedures lead to similar values of the renormalized charge. In the case of 



salt-free suspensions, the effective charge is found to saturate at [22 



Z'rt ^ , (125) 

where x is an approximately linearly increasing function of colloidal concentration for 
suspensions with volume fraction larger than 1%. For suspensions with colloidal volume 



fraction between 1% and 10% the value of x varies from around 9 to 15 ||112| , |113 |. 



5.5. Colloidal Fluid 

In this section we will apply the insights gained from the study of one and two component 
plasmas to the exploration of stability of charged colloidal suspensions against a gas- 
liquid phase separations. We note that the large size asymmetry between colloids and 
counterions leads to very different equilibration time scales. On the time scale of polyion 
motion, the counterions are always equilibrated. This suggests that the calculation of 
free energy should be done in two stages ||114|| . First, we shall trace out the counterion 



degrees of freedom, leading to effective many-body interactions between the colloidal 
particles. Then we will use these effective interactions to calculate the colloid-colloid 
contribution to the total free energy. The procedure is similar to the one used in 



McMillan-Mayer theory of solutions [p.l5|| . 



We shall ffist calculate the contribution to the total free energy arising from the 



polyion-counterion interactions Consider a suspension in thermal equilibrium. 

While the colloidal particles are more or less uniformly distributed throughout the 
solution, the positions of counterions are strongly correlated with the positions of 
polyions. As a leading order approximation we can, therefore, take the polyion-polyion 
correlation function to be 

9pp = 1 (126) 



32 



Electrostatic correlations: from Plasma to Biology 



and the polyion-counterion correlation function to be 

V = e-^^'^(^) . (127) 

Choosing the coordinate system in such a way that it is centered on top of one of the 
colloidal particles, the electrostatic potential at distance r < a satisfies the Laplace 
equation, while for distances r > a it satisfies the Poisson equation Eq. (^. Based on 
Eqs. (T^) and ( |127| ) the charge density in the region r > a can be approximated by, 

p,(r) = -Zgpp + gp,e-^^*W . (128) 

In the spirit of the Debye-Hiickel theory we shall linearize the exponential [^, H. The 
distribution of charge around the colloid reduces to 

P,{r) = -f3q'pc<P{r) • (129) 

For r > a the electrostatic potential, therefore, satisfies the Helmholtz equation (|^) with 
K given by Eq. (|113|) . The solution to this equation is 

0>W = , (130) 

er 

while the solution to the Laplace equation for r < a is 

'^<W = ■ (131) 

ea(l + na) 

The electrostatic energy due to polyion-counterion interaction is 

Up = \j dh[p,{v) + gp(r)]0(r), (132) 

where Pg(r) is the charge density of counterions given by Eq. ( |129| ), and qp(j) is the 
charge density of a polyion, 

gp(r) = -^5(|r|-a). (133) 
Performing the integration we find 



2„2 

(134) 



2e(l + Ka 



1 



a 2(1 + Ka 

The electrostatic free energy of a polyion inside the suspension is obtained using the 
Debye charging process ||116||, 



Tp = / dX = . 135 

Jo A 2ea(l + Ka) 

Note that this free energy is the sum of the polyion self energy 

r-/f = ^, (136) 



and the solvation energy 

csolv 

^ ~ 2ea(l + Ka 



jrsoi. ^ _ zyKa ^ ^^g^^ 
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which the polyion gains from being inside the "ionic sea" . The electrostatic free energy 
due to interaction between all the polyions and counterions is 

F-=-^!^!^. (138) 
2ea(l + Ka) ^ ^ 

We have effectively integrated out the counterion degrees of freedom. This, however, 
leaves us with the effective many-body potentials of interaction between the colloidal 
particles. For dilute suspensions, the pairwise interaction potential should be the 
dominant one. The two-body interaction potential can be obtained from the solution 
of Helmholtz equation for two colloidal particles ||117| , |118|| . At large separations 
this leads directly to the DLVO interaction potential, Eq. (|115|) . This potential has 
been extensively tested experimentally and found to work very well for bulk colloidal 
suspensions ||119|| . Since the DLVO potential is short ranged, the contribution to the 
total free energy arising from the colloid-colloid interaction can be calculated in the 
spirit of the traditional van der Waals theory, through the second virial term. A more 
sophisticated calculation of the colloid-colloid free energy relies on the Gibbs-Bogoliubov 
variational bound, 

ppp <Fo + {Vdlvo)o , (139) 

where the reference system is taken to be the fluid of hard spheres, whose diameter plays 
the role of a variational parameter. The free energy resulting from the polyion-polyion 
interaction, F^p, can be approximated by the lowest variational bound of Eq. (|139|) . 
The calculation is somewhat involved, so we refer the interested reader to the original 



papers M 121, 20, 122 



The entropic mixing free energy of colloids and their counterions is simply that of 
an ideal gas, 

= ZNp[\n{Zpj,Al) - 1] + Nj,[\n{pj,Al) - 1] , (140) 

where Ac and Ap are the de Broglie thermal wavelengths of counterions and polyions, 
respectively. 

The total free energy of colloidal suspension is the free energy needed to solvate 
colloids in the sea of other polyions and counterions F^'^ + F^^, and the free energy of 
mixing 

p ^ ppc _^ ppp pent _ {I'^l) 

The osmotic pressure is 
dF 

P = . (142) 

dV N, ^ ^ 

It is found that for suspensions with 

C = ±^ > 15.2 (143) 
a 

the pressure is not a convex function of the colloidal concentration, implying existence 
of a thermodynamic instability. At criticality the colloidal volume fraction is around 
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1%. The crucial question is whether this result is reliable? In order to calculate the 
electrostatic free energies, we were forced to linearize the Boltzmann factor. While this 
is a reasonable approximation at large separations away from the polyions, linearization 
is clearly invalid in the vicinity of colloidal surface. There, the strong electrostatic 
interactions result in an accumulation of counterions and the effective polyion charge 
renormalization. Therefore, the linear theory can be used only if the bare colloidal charge 
is replaced by the effective renormalized charge, Z ^eff, in all the expressions. It 
was found, however, that the bare charge does not increase without limit but saturates 
at the value given by the Eq. (|125|) . Substituting Z ^e//, into the definition of C 



Eq. (|143D , we see that C < 15 for all the values of the bare charge Z in the critical 
region. The critical threshold, therefore, can not be reached, meaning that a deionized 
aqueous suspensions with monovalent counterions is stable against a liquid-gas phase 
separation for all colloidal charges and sizes. This conclusion has also been confirmed 



by more detailed calculations and simulations p3|, |123| , p^ , |124 . 

The result that the non-linear terms omitted within the Debye-Hiickel 
approximation stabilize a deionized colloidal suspensions against a liquid-vapor phase 



separation has also been obtained by von Griinberg et al. ||125| , |126| , [127|| and Tamashiro 



and Schiessel [|128|| based on the analysis of the full non-linear Poisson-Boltzmann 
equation inside a Wigner-Seitz cell. The numerical integration of the non-linear PB 
shows that the osmotic pressure is a monotonically increasing function of colloidal 
concentration. This means that at the level of WS approximation suspension is 
thermodynamically stable. Von Griinberg et al. and Tamashiro and Schiessel, however, 
demonstrate that the linearized PB equation leads to the negative compressibility and 
the osmotic pressures for highly charged colloidal particles. This erroneously suggests 
presence of a thermodynamic instability. Clearly the instability is an artifact of the 
linearization. Furthermore, our calculations show that any linear theory, which does not 
take into account the colloidal charge renormalization, is likely to lead to an incorrect 



prediction of a liquid- vapor phase separation pO|, p.29|| in deionized aqueous suspensions 
with monovalent counterions. 

It is curious that the "linear" correlations between the colloids and the counterions, 
responsible for the screening of electrostatic interactions, are also the ones driving 
the suspension towards the phase separation. On the other hand, the "non-linear" 
correlations responsible for the counterion condensation and the colloidal charge 
renormalization, stabilize the suspension against a phase transition. 

6. Polyelectrolyte solutions 

Polyelectrolytes are polymers with ionizable groups which have tendency to dissociate 
in polar solvents | |130| , |131|| . The good water solubility of polyelectrolytes is due to large 



favorable gain in the solvation free energy resulting from hydration of charged monomers 
and counterions. Unlike polyampholytes ||132| , |133| , |134| , |135|| , whose monomers can either 



be cationic (positive) or anionic (negative), all charged monomers of a polyelectrolyte 
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carry charge of the same sign. Depending on the sign of this charge, polyions are 
either cationic or anionic. Over the last few decades polyelectrolytes have found many 
industrial applications ranging from water treatment and oil recovery to detergents 
and superabsorbants. The biological importance of polyelectrolytes, however, has been 
realized much earlier. After all the most important biomolecule, DNA, is an anionic 
polyelectrolyte ||106 |. 

Unlike the simple polymeric fluids, thermodynamics of which is fairly well 
understood, polyelectrolyte solutions still remain to large extent enigmatic. The 
difficulty in studying polyelectrolytes resides in the combination of polyion 
flexibility [|136|, H^, HH, |130|, |13|, |4|, [l4ll, the long-ranged nature of the Coulomb 
force, as well as a large charge and size asymmetry between the polyions and the 
counterions. There are, however, some polyelectrolytes whose polyions are rigid 
molecules. This allows to bypass the complications associated with the statistics of 
polyion conformations. We have already come across this kind of systems in our 
exploration of charged colloidal suspensions. Many biologically relevant polyelectrolytes 
are also fairly rigid. Persistence length (the distance over which polymer can be 
considered to be rodlike) for double stranded DNA is on the order of 500 A. Actin 
filaments, which are the building blocks of a cytoskeleton, have persistence length even 
larger, on the order of microns. This should be compared with the Debye length at 
physiological concentrations of 150 mM of NaCl, which is about 8 A. Clearly under 
these conditions the flexibility of the DNA or the actin fllaments can be considered an 
irrelevant perturbation. 

The thermodynamics of rodlike polyelectrolytes can be explored using the same 
theoretical tools used for spherical colloidal suspensions. Rodhke molecules can undergo 
nematic and smectic phase transitions. For ordered periodic structures, the Wigner-Seitz 



cell formalism can be employed to obtain most of the relevant thermodynamics ||142 



143|| . At low volume fraction, when a polyelectrolyte solution is disordered, this strategy 



is no longer valid and a different methodology must be used. This can be constructed 



along the same lines taken for colloidal suspensions [144, 145 1. The fundamental role 



of electrostatic correlations between the polyions and the counterions appears as Debye 
screening of polyion-polyion interactions and the renormalization of a polyion charge. 
In polyelectrolyte literature the polyion-counterion association leading to polyion charge 
renormalization is known as the Manning condensation ||104| , |106| , |105|| . Here we will 
show that the Manning condensation is very similar to the charge renormalization found 
in colloidal suspensions. 

First, we shall briefly review Manning's original argument ||104|| . Manning was 



interested in deriving the limiting (low density) laws for polyelectrolyte solutions, similar 
to the ones found by Debye and Hiickel for simple electrolytes, see Eq. (0). The salient 
feature of the Debye-Hiickel limiting laws is that they do not depend on specifics of 
electrolyte, i.e. size or hydration. For example, the osmotic pressure at low ionic 
concentrations is found to be a function only of the ionic charge, temperature, and 
concentration. The question then arises if such a limiting law is also possible for 
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polyelectrolyte solutions. The answer to this question is far from obvious. The strong 
electrostatic interaction between the polyions and counterions favors accumulation of 
counterions in the vicinity of polyions. It is, therefore, possible that even at very large 
dilutions the physics of a polyelectrolyte solution remains that of a strongly interacting 
system for which no limiting law should be anticipated ||106 . 



6.1. Manning condensation 

Consider a simple model of a polyelectrolyte solution. The rodlike polyions of 
concentration pp, idealized as rigid cylinders of length L and radius a, carrying Z ionized 
groups — each of charge —q uniformly spaced along the major axis of the cylinder — 
inside a uniform dielectric solvent of constant e. The counterions of concentration 
Pc = Zpp will be treated as point particles of charge q. For simplicity we will restrict 
our attention to the situation in which there is no added salt. 

In the low density limit we can neglect the discreteness of the polyion charge 
distribution and assume a uniform line-charge density, 

A„^-f.-f, (144) 

where b is the separation between the successive charged monomers along the 
polyelectrolyte chain. The bare interaction potential between a long charged cylinder 
and a counterion is, 

,145) 

where tq is the arbitrarily chosen point of zero potential. A polyion-counterion two-body 
partition function is 

J a 1 S 

where R is the cutoff distance at which a counterion cas still be considered to be bound 
to the polyion. The Manning parameter is defined as, 

« = lf^ = ^. (147) 



The integral in Eq. ( |146|) remains finite for all values of ^. Manning noticed, 
however, that if the limiting laws exist, the thermodynamic functions should be 
independent of the polyion diameter. However, if a = the integral in Eq. ( |146| ) diverges 
as ^ —>■ 1~. Manning interpreted this divergence as an indication of the counterion 
condensation. For values of ^ > 1 he supposed that n counterions condense onto the 
polyion, reducing proportionately its effective charge density from Aq to 

An = • (148) 

To find the number of condensed counterions n, Manning postulated that for ^ > 1 the 
effective reduced line charge density 
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saturates at one, 

^eff = 1 . (150) 

If the "renormalized" ^e// is used in Eq. ( p.46| ), instead of the "bare" ^ when ^ > 1, the 



polyion-counterion partition function remains finite. Eq. ( |149| ) and Eq. ( [150| ) determine 
the number of condensed counterions to be 

n* = Z (l-l] for i>l 



= for i<l (151) 

Once n* is determined, the rest of the thermodynamic functions can be calculated 
quite easily. The nice thing about Manning's argument is that it is so simple. On the 
other hand it contains some points which might leave a more mathematically inclined 
reader quite disturbed. Manning relied on the existence of the limiting laws to establish 
the limiting law in the first place! This is a circular logic which is not always guaranteed 
to work. It is interesting to apply the same argument to the case of a two-dimensional 
plasma of particles interacting by a logarithmic potential [ p,46|| . Section ^. Consider 
an anion-cation two body partition function Eq. (|39|). Following Manning, lets look 
at the limit a — > 0. In this case we find that the integral in Eq. (^) diverges for the 
temperatures 

T<r..^. (152) 

Thus, we might incorrectly conclude that the metal-insulator transition also happens 
at Td- In reality, we know that the Kosterlitz-Thouless transition occurs at half this 
temperature, 

r = r...jlj. (153) 

For any finite value of the particle diameter, all the thermodynamic functions are 
analytic at Td and the singularities only appear at T = Tkt- Therefore, it seems 
far from obvious if Manning's argument is valid for real polyelectrolytes with monomers 
and counterions of finite diameter. To explore this in more detail we can appeal to 
the Debye-Hiickel-Bjerrum theory, previously constructed for low dielectric electrolytes. 
Section O. 



6.2. Counterion Association 

We will restrict our attention to the salt-free infinite dilution limit. The calculations, 
however, can be extended to solutions of finite polyelectrolyte concentration as well as 
to polyelectrolytes with salt or even amphiphiles ||145| , |147|| . As in the case of colloidal 
fiuids, we would like to trace out the degrees of freedom associated with the counterions. 
This leads to the effective many-body interactions between the polyions. In the limit of 
infinite dilution the contribution from these interactions to the total free energy can be 
neglected. 
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Consider a dilute polyelectrolyte solution in thermal equilibrium. If we take 
a snapshot, we will see polyions distributed more or less uniformly throughout the 
solution, with no specific orientation. The counterions, on the other hand, will be 
clustered in the vicinity of polyions. This picture suggests that there are only weak 
positional correlations between the polyions and strong positional correlations between 
the polyions and the counterions ||148|| . The average distribution of charge around a 



polyion can, therefore, be approximated by the Eq. (|128 



Suppose we choose a coordinate system so that it is centered on one of the 
polyions, with the z-axis along the polyion's major axis. The electrostatic potential 
for r < a satisfies the Laplace equation, while for r > a it satisfies the Poisson 
equation with the charge distribution approximated by the Eq. ( |128| ). As in the case of 



colloids, we would like to linearize the exponential Boltzmann factor. This, however, is 
prohibited by the fact that electrostatic interactions are very strong in the vicinity of the 
polyion surface. Thus, in order to linearize the exponential, the short distance polyion- 
counterion correlations must be explicitly taken into account. In our earlier study of 
colloidal suspensions this was done by introducing an effective renormalized charge. Here 
we will take a somewhat different approach. The main consequence of short distance 
electrostatic correlations is the polyion- counterion association. This is very similar to 
the concept of Bjerrum association, which has proven so successful for symmetric 1:1 
electrolytes. A polyelectrolyte solution can be thought of as being composed of free 
unassociated counterions of density pf, as well as of clusters of density p„, consisting of 
one polyion and some number < n < Z of associated counterions. The polyions with 
no associated counterions are treated as 0-clusters. The particle conservation requires 
that 

z 

^Pn = Pp. (154) 



n=0 



and 



Pf + ^^Pn = Zpp . (155) 
n=0 

The distribution of cluster sizes {pn} can be determined from the equilibrium condition 
that the Helmholtz free energy be minimum. Once the clusters are explicitly introduced 
into the theory, the exponential factor in Eq. ( |128|) can be linearized, since at large 
distances (3q(f){r) < 1 and at short distances the "non-linearities" are accounted for 
through the cluster formation. After linearization, the Poisson-Boltzmann equation 
reduces to the Helmholtz equation with k given by. 



The linear equation can be easily solved yielding the electrostatic potential outside a 
n-cluster 

2A„ KQ{Kr) 

0(r) = , . , (157) 

e KaKi[Ka) 
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where A„ is the linear cluster charge density of a n-cluster Eq. ( |148| ), and K^{x) is the 
modified Bessel function of order u. For distances r < a the electrostatic potential is 
found to be 

-K'o(Ka) 



0(r) = ^ln(r/a) + " 



(158) 



e ■ ■ e KaKi[Ka) 
The electrostatic energy due to n-cluster-counterion interaction can be obtained from 
Eq. (|132|) and the electrostatic free energy follows from the Debye charging process as 
in Eq. (|135|) . In the limit of infinite dilution, the electrostatic free energy of a n-cluster 
is 



n] 



27i 



+ ln 



(f) 



(159) 



The electrostatic free energy density / = F/V due to all cluster-counterion interactions 
is 



Pf 



pc 



-1^^{Z -nfpr, 



n=0 



27J 



— + !„(-) 



(160) 



The cluster-cluster and the counterion-counterion contributions are of higher order in 
density and in the limit of infinite dilution can be neglected. The only contributions 
which must still be taken into account are the entropic free energy of mixing and the 
free energy necessary to construct the isolated clusters. Both of these can be concisely 
written as 

'.3(n+l), 

2^p- 



cent—cl 



In 



Pn 



Cn 



- 1 



(161) 



where the de Broglie thermal wavelength of a n-clusters is 
A. 



and rrin is the cluster geometric mean mass, 



rrir, 



The internal partition function of a n-cluster is 



Cn I 

nl 



i=l 



-I3U 



(162) 



(163) 



(164) 



were U is the usual Coulomb potential. A suitable cutoff must be chosen in order to 
define what constitutes a cluster. Evaluation of the integral in Eq. ( p.64| ) represents a 
formidable task. Fortunately, as we shall see, for polyelectrolyte solutions at infinite 
dilution specific knowledge of Cn proves unnecessary. The total free energy density is 
then 



fiipn}) = + r 



ent—cl 



(165) 
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To find the equilibrium cluster distribution, this free energy must be minimized subject 
to constraints of particle conservation, Eqs. ( p.54| ) and ( |155| ). The minimization is 
equivalent to the law of mass action, 

IJ.n = fJ'O + niJf , (166) 

where the chemical potential of n-clusters is 
df 

f^n= , (167) 

OPn 

and the chemical potential of free ions is 

^'=§-r (168) 

Substitution of the total free energy into Eq. ( |166|) leads to the n-cluster distribution 

Pn = CnPoP/e^^''°^+"^^^^-'^"^^ , (169) 

which is actually a set of Z coupled algebraic equations. At large dilutions, the excess 
chemical potential of n-clusters is 

and of free ions is 

m' = -it^-^^£^ + 0(,p,). (171) 

The internal partition function of a ra-cluster is independent of density. Recalling 
that K ~ a/P/, Eq. (|169|) simplifies to 

Pn ~ Popf^ , (172) 
where the exponent is 

gin) = j^n^ - + n , (173) 
In the limit of infinite dilution, — 0, only the cluster of size 

nM = z(l-l] , (174) 



OiPf) , (170) 



which minimizes g{n) survives. In this limit the cluster size distribution takes a 
particularly simple form, 

Pn = PpSnnM ■ (175) 

This is precisely the cluster size distribution postulated by Manning based on his 
heuristic argument. The osmotic pressure of a polyelectrolyte solution can be obtained 
through the Legendre transform of the negative free energy density | |146| |, 

P = -fipfApn}) + PfPf + ^Pnpn , (176) 
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PP={l-^]Zpp for e<l- (177) 



which leads directly to the Manning limiting law for pressure [ [104 

1 

PP=^ for e>l. (178) 

The discontinuity in the slope of pressure as a function of temperature has 
provoked a lot of speculation that the Manning condensation is a real thermodynamic 
phase transition. Kholodenko and Beyerlein went even so far as to identify Manning 
condensation with the Kosterlitz-Thouless transition ||149|] . That this is incorrect 



follows already from our discussion in Section |6.1| . For the two dimensional plasma 
with logarithmic interactions, the Kosterlitz-Thouless transition occurs at half the 
equivalent Manning temperature. The two phenomena, therefore, have nothing in 



common ||146|| . Furthermore, while the Kosterlitz-Thouless is a real thermodynamic 



phase transition characterized by the diverging Debye length, for rodlike polyelectrolytes 
the discontinuity in the slope of pressure appears only in the double limit Pp — > 0, 
L —>■ oo. If the order of limits is interchanged, there is no singularity and no counterion 
condensation. In addition directly at the condensation threshold ^ = 1, the Debye 
length remains finite [ |145| . 



For polyelectrolyte solutions containing polyions of finite size and at non-zero 
concentration, there is also a polyion-counterion association ||145|| . In this case, however, 
the distribution of cluster sizes is no longer a delta-function, but rather a bell-shaped 
curve centered on the value n*. We find that n* depends on the concentration of 
polyelectrolyte and is somewhat larger than the limiting Manning value um- The 
pressure remains an analytic function of ^, showing that for real polyelectrolyte solutions, 
the counterion condensation is actually a crossover phenomenon, very similar to the 
micellar formation in amphiphilic systems [|150 . 



7. Multivalent counterions 

Up to now we have been concentrating our attention on aqueous solutions with 
monovalent counterions. It was already mentioned that in this case the correlations 
between the condensed counterions can be neglected. To understand why, let us 
compare the characteristic electrostatic energy of a counterion-counterion interaction 
to the characteristic thermal energy ksT, 

- (179) 



edhsT 

where a is the counterion valence and d is the average separation between the n 
condensed counterions on the surface of a colloidal particle of radius a. Since nvrrf^ = 
47ra^, 

2a , , 

(180) 
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and the coupling strength becomes, 

r = . (181) 

Now, lets consider highly charged latex particles with Z = 7000 and a = 1000 A, in 
water at room temperature. From Eq. ( 125 ), taking x = 15, Z^j^j- = 2100, which means 



that 4900 monovalent {a = 1) counterions are condensed onto the particle. The coupling 
strength of the counterion-counterion interaction is then F ^ 0.25, which clearly shows 
that the electrostatic interactions between the condensed counterions are very weak. We 
can make this observation even more general. The high surface charge concentration a^n 
encountered in nature is on the order of one elementary charge per 100 A2. Lets suppose 
that suspension consists of highly charged colloidal particles with surface charge density 
am- Clearly this means that there will be a lot of counterion condensation. For a salt- 
free colloidal suspension containing multivalent counterions, the number of condensed 
counterions will be approximately 

n* ^- . (182) 

a 

The radius of a colloidal particle can be expressed as 

(183) 



47rcr- 



m 



Substituting Eqs. ( |182| ) and (|183|) into Eq. (|181|) we find that the maximum counterion- 



counterion coupling strength is, 

3 

Tmax ~ a^Xsy/T^ ■ (184) 

For monovalent counterions Tmax ~ 1-3, for divalent counterions Tmax ~ 3.6, and for 
trivalent counterions Tmax ~ 6.8. Although Tmax is an overestimate, it clearly shows that 
for highly charged colloidal particles, correlations between the condensed multivalent 
counterions cannot be ignored. 

7. 1 . Overcharging 

One consequence of strong electrostatic correlations is the phenomenon known as the 
"overcharging" |p.51| , |152| , |153| , |154| , |155| , |156| , |157| , |158| , p.59| , |160|| . Overcharging occurs 



as the result of highly favorable gain in electrostatic free energy due to strong positional 
correlations between the condensed counterions. 

To understand better how the overcharging of colloidal particles comes about let 
us consider a simple case of one colloidal particle with a uniform surface charge —Zq 
and radius a, at zero temperature ||161| , |162|| . The question that we would like to answer 



is how many a-valent counterions should be placed on top of the colloidal particle in 
order to minimize the electrostatic energy of the resultant polyion-counterion complex? 
Naively we might suppose that the number of condensed counterions should be such 
as to neutralize completely the colloidal charge. This, indeed, would be the case if 
the charge of counterions was uniformly smeared over the surface of colloid. In reality, 
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the counterions are discrete entities and can gain favorable energy by maximizing their 
separation from one another. Lets calculate the electrostatic energy of the polyion- 
counterion complex, 

2ea ea 

The first term is the self energy of a polyion, the second term is the electrostatic energy 
of interaction between the polyion and n condensed a-ions, and the last term is the 
electrostatic energy of repulsion between the condensed counterions. Now, consider a 
one component plasma of n a-ions on the surface of a sphere of radius a but with a 
uniform neutralizing background charge —anq. The electrostatic energy of this OCP 
can be expressed as the sum of contributions arising from the counterion-counterion 
interaction, counterion-background interaction, and the self energy of the background, 

^ocP^^„._^^^, (186) 

ea lea 

Substituting this expression into Eq. ( |185|) , the electrostatic energy of the polyion- 
counterion complex simplifies to 

For low temperatures, the condensed counterions try to maximize their separation from 
one another. In the planar geometry the ground state corresponds to a triangular 
Wigner crystal. A similar arrangement of counterions will also be found on the surface 
of a spherical colloidal particle, up to some topological defects. The electrostatic energy 
of a planar OCP has been discussed in Section For a spherical OCP the electrostatic 
energy at zero temperature is 

F^CP = _M^!^^ . (188) 
2ea 

where M is the Madelung constant. Because of the topological difference between the 
plane and the surface of a sphere, we expect that the Madelung constant will not be 
exactly the same in the two cases. The difference, however, should not be very large as 
was confirmed in recent Monte Carlo simulations | [L62|| . For concreteness, we shall use 
M = 1.106, the value of the planar OCP. 

The effective charge of the polyion-counterion complex in units of —q is 

Z^fj = Z -an* . (189) 

where n* is the number of condensed a-ions which minimize the electrostatic energy, 

^ = . (190) 

an n* 

The effective charge is found to be 

Z - 1 + ..... 
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where 7 is, 

We see that the effective charge of the polyion-counterion complex is inverted compared 
to the bare charge Z of the colloidal particle, so that the complex is overcharged. 
For highly charged colloids, the effective charge scales as the square root of the bare 
charge [|163], |16|], 



Zeff ^ . (193) 

The analysis above was conducted for one colloidal particle at zero temperature. For 
solutions at finite concentration and temperature we face the same difficulties already 
encountered in our earlier discussion of charge renormalization in colloidal suspensions 
and polyelectrolyte solutions. In fact the problems are even more severe, since for 
multivalent counterions the Poisson-Boltzmann equation fails completely. A logical 
step is to appeal to the density functional theory. The difficulty with the DFT lies in 
constructing a suitable density functional which can take into account the electrostatic 
correlations. For Z : a suspensions without salt, such functional was presented in 
Section |5]^. It was found using the DFT and the Wigner-Seitz cell formalism pl[] , 
that the effective charge Z^ff, as the function of the bare charge Z, after reaching the 
maximum decreases, asymptotically going to zero as Z ^ 00. This behavior is in striking 
contrast to the saturation of the effective charge predicted by the Poisson-Boltzmann 
theory [^. Unfortunately it is difficult to construct a suitable density functional for 
suspensions containing salt. One alternative is to use the integral equations. The 
numerical complexity of these theories, however, tends to obscure the essential physics of 
the problem. Below we shall present a simple phenomenological model of overcharging. 
Our goal is not the quantitative accuracy, but rather the physical insight into the 
mechanisms leading to the overcharging in the polyelectrolyte solutions. 



7.2. Overcharging in electrolyte solutions 

Consider a dilute colloidal suspension containing a monovalent salt at concentration 
c, and an a-valent salt at concentration c^. In aqueous solution the monovalent salt 
dissociates producing 1 : 1 electrolyte, while the a-valent salt dissociates into a : 1 
electrolyte. The inverse Debye length is 

= K = ^StiXbI , (194) 

where the ionic strength is 

I = — [a^Ca + aCa + 2c) . (195) 

For simplicity we shall restrict our attention to suspensions with monovalent salt 
near physiological concentrations — 150mM of NaCl. Under these conditions the 
Debye length is around 8 A, and the interactions between the colloidal particles can 
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be completely neglected. Furthermore, since the electrostatic attraction between the 
highly charged a-ions and colloids is so much stronger than the interaction between 
the monovalent counterions and colloids, the effective charge of the polyion-counterion 
complex is completely determined by the number of condensed a-ions. 

We will define the counterions as free (not-associated) if they are farther than some 
distance S from the colloidal surface. The "agglomerate" is then defined as the polyion 
with its 5-sheath of surrounding counterions. We note, however, that not all of the 
a-ions in the agglomerate are actually associated with the polyion. The reason for this 
is that many of these ions have sufficient kinetic (thermal) energy to leave the vicinity 
of colloidal particle. Only the counterions which have the total (kinetic plus potential) 
energy less than zero can be considered bound to the polyion. Unfortunately, it is not 
easy to come up with a practical implementation of this energetic criterion, except within 
a Molecular Dynamics Simulation. On the other hand, a simple geometric criterion 
based on the polyion-counterion separation is easily implemented. Care, however, must 
be taken not to count all of the a-ions inside the agglomerate as belonging to the 
polyion-counterion complex. Bellow we shall see how this can be accomplished. 

We shall take 6 to be on the order of few angstroms, corresponding to the radius 
of a hydrated ion, 5 ~ 2 A. Since the agglomerate is in contact with the bulk, its size is 
determined by the minimum of the grand potential function 

Q{n) = F{n) - n/i„ , (196) 

where F{n) is the Helmholtz free energy of the agglomerate and /Iq, is the chemical 
potential of a-ions inside the sheath. In thermal equilibrium, the chemical potential 
of a-ions inside the sheath equals to the chemical potential of a-ions in the bulk of 
the suspension. For low bulk concentrations, /i^ can be approximated by the chemical 
potential of an ideal gas, 

/i„ = ln(c«A^) . (197) 

The Helmholtz free energy of an agglomerate is then 

Fn = En + + F^* . (198) 

The electrostatic free energy of an isolated polyion-counterion agglomerate E^ is given 
by the Eq. ( |187|) , the solvation energy that the agglomerate gains when placed in an 
ionic environment is, 

2ea(l + na) 

see Eq. ( |137| ), and the entropic energy of counterions inside the sheath is 

Fr* = fcBT[nln(p„A=')-n] . (200) 

where, 
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For high valence counterions (strong-coupling limit) the free energy of the OCP can 
be approximated by that of a Wigner crystal, see Section H. If more accuracy is 



needed, one can use the extrapolation formulas based on Monte Carlo simulations |80 
Here, however, we shall content ourselves with the simplest approximation. The grand 
potential of an agglomerate is, 

n) = 4^^^^ - M^^!^ + n HpJc.) - n . (202) 
2a(l + Ka) 2a 

The number of a-ions n* inside an agglomerate is determined from the minimum of the 

grand potential, 

an(z, n) 

dn 

From the previous discussion recall that it is incorrect to associate the effective charge 
of the polyion-counterion complex with the value of n*, i.e. Zeff Z — an*. The reason 
for this is that not all of the a-ions inside the 5-sheath are actually bound to the polyion. 
The real number of condensed counterions is n* — n*, where the overestimate ra* can be 
obtained by considering the number of a-ion within the distance 6 from the surface of 
a "neutral" polyion, Z = 



. (203) 



dn{0,n) 



. (204) 



dn 

The effective charge of the polyion-a-ion complex is then 

Zeff = Z - a{n* - <) . (205) 

If a small electric field is applied to the suspension, it is the Z^ff which will determines 
the electrophoretic mobility of colloidal particles ||160| , |110| , |164| , 

In Figs. § and |] we present the effective colloidal charge as a function of 
concentration of divalent and trivalent counterions, for suspension containing colloidal 
particles with Z = 4000 and a = 1000 A. It is curious to note the appearance of a 
minimum in the effective charge as a function of trivalent ion concentration. In Fig. ^ 
we present a plot of the effective charge as a function of bare charge for colloids with 
a = 1000 A , in a suspension containing monovalent salt at concentration c = 0.15M 
and trivalent counterions at C3 = O.OIM. We note that unlike suspensions containing 
monovalent counterions, the effective charge does not saturate, instead it reaches the 
maximum value and then falls off sharply. For colloids with Z ^ 11500, the bare 
colloidal charge is completely neutralized by the a-ion condensation. In Fig. |^ we show 
the dependence of Zeff on the amount of monovalent salt. For small concentrations of 
a-ions, the effective charge of the complex is found to increase with the concentration 
of monovalent salt, asymptotically approaching the bare value. However, when the 
concentration of multivalent ions reaches the critical value, there is a qualitative change 
of behavior. At this point the effective charge is no longer a monotonically increasing 
function of the monovalent salt concentration. Instead after reaching the maximum, 
Zeff begins to decline, eventually going through the isolectric point {Zeff = 0) and 
charge inversion. Fig. |^. 
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Figure 4. The effective (renormaiized) cfiarge of colloidal particles oi Z = 4000, a = 
1000 A inside a suspension containing monovalent salt at physiological concentration 
of 0.15M, as a function of concentration of divalent counterions. 
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Figure 5. The effective (renormaiized) charge of colloidal particles oi Z = 4000, a = 
1000 A inside a suspension containing monovalent salt at physiological concentration 
of 0.15M, as a function of concentration of trivalent counterions. 
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Figure 6. The effective (renormalized) cliarge of a colloidal particle with a = 1000 A 
as a function of the bare charge Z, inside a suspension containing monovalent salt at 
physiological concentration of c = 0.15M and the trivalent ions at C3 = O.OIM. 




Figure 7. The effective (renormalized) charge of colloidal particles with Z = 4000 
and a = 1000 A inside suspension containing trivalent counterions with C3 = O.OIM, 
as a function of concentration of monovalent salt c. 
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Figure 8. The effective (renormalized) charge of cohoidal particles with Z = 4000 and 
a — 1000 A inside the suspension containing trivalent counterions with C3 = 0.0135M, 
as a function of concentration of monovalent salt c. 



We have considered only the simplest form of overcharging involving colloids and 
multivalent microions. There is a number of interesting variations. Recent developments 
in the field of gene therapy require construction of safe and efficient trans-cellular gene 
delivery systems ||165| , p.66| , |167]| . Both the DNA and the phospholipids, which are the 
main constitutive component of a cellular membrane, are negatively charged. There is, 
therefore, a strong electrostatic repulsion between the DNA and the cellular membrane. 
This repulsion inhibits the transfection of naked DNA into the cells. Furthermore, in 



vivo the unprotected DNA is rapidly degraded by the nucleases present in plasma [|168 

Much of the effort of gene therapy has been concentrated on viral transfection. 
A retrovirus (virus which incorporates its genetic material into the host genome) or 
adenovirus (which does not) has its genetic material removed and substituted by the 
gene that needs to be replicated. The modified viruses are then made to infect the 
cells, thus effectuating transfer of the genetic material. There are, unfortunately, a 
number of serious complications involved with this procedure. These range from a strong 
immunological response of an organism against the infecting virus, to potentially deadly 
consequences arising from the recombinant viral structures ||168|| . All these factors have 
lead to attempts to develop non-viral transfection methods. One of the promising 
approaches relies on formation of the DNA-cationic lipid complexes, or lipoplexes for 
short. The hydrophobic interaction between the lipid tails, in addition to the cationic 
charge of their head groups, favors their agglomeration in the vicinity of polyions. For 
sufficiently long hydrocarbon chains, the gain in the hydrophobic energy resulting from 
the lipid condensation onto the DNA is sufficient to lead to the charge reversal of 



a lipoplex ||169||. This happens even though the charged head group of the lipid is 
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monovalent. The overcharging in this case is the consequence of hydrophobicity of lipid 
molecules and is not the result of electrostatic correlations. 

A different method for gene transfection uses multivalent counterions such as Ca^^ . 
The efficiency of Ca"*""*" as a transfection agent might be due to its ability to neutralize 
or even invert the helical DNA charge without making it collapse. In this respect it is 
quite different from the transition metal ions such as Mn"^"*" and Crf^"*" which neutralize 
the DNA charge, but also lead to its condensation [|170|| . Formation of a neutral or 



overcharged DNA-Ca++ complex allows the DNA to come into a close proximity of the 
cellular membrane. Transfection might then be able to proceed through a reptation-like 
motion of the DNA through a cellular pore. 

A very curious form of overcharging is found to occur in cellular compaction of the 
DNA. Virtually all of the DNA in nuclei of eukaryotic cells exists as a highly organized 
nucleoprotein fiber called chromatin. At the lowest level of the chromatin hierarchical 
structure is a nucleosome. Nucleosome resembles a thread wound around a cylindrical 
spool. The thread is the DNA molecule, while the spool is an octameric protein 
composed of eight smaller proteins called histones. Each octamer has 220 basic residues 



of which 148 are on the surface of the protein and are exposed to the solvent ||171|| . The 
rest are inside the protein core and are unlikely to be ionized. The maximum charge of 
an octamer is, therefore, Q^^^ = +220e, but is likely to be significantly lower than this 
value under the normal physiological conditions. Each octamer is encircled by 1.8 turns 
(147 base pairs) of the DNA thread, carrying a net charge of — 294e. A nucleosome is, 
therefore, strongly overcharged by the associated DNA [|157|, |172|, |173|, |174 . 



It is interesting to speculate how the nature is using this overcharging in the overall 
organization of the chromatin. The degree of compaction achieved inside the nucleus is 
quite astonishing. The total length of the DNA in the nucleus of a human cell is about 
3.3 billion base pairs. If extended it would be more than Im long. Yet, it is compacted 
into a nucleus of diameter of lO/im! Furthermore, all this compaction is done in such a 
way that the DNA is easily transcribed and replicated. This is, indeed, an astonishing 
feat of chemical engineering! 

8. Like-charge attraction 
8.1. Confined suspensions 

One of the most curious phenomena which has produced much healthy debate in the 
condensed matter community is the appearance, under some circumstances, of attraction 
between like-charged macromolecules. The first observation that something might be 
missing in the traditional DLVO theory of colloidal stability came from the experimental 
observations of Ise et al. of void structures inside the highly deionized suspensions [p.02| . 
The finding of voids and clusters has lead Sogami and Ise to propose a modification of 
the DLVO pair potential. In its stead they suggested a new potential obtained from 
the considerations of the Gibbs free energy of suspension. The Sogami-Ise potential 
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has created a lot of stimulating controversy, which is still in progress. The problem 
has been re-analyzed by Overbeek ||103|| , who has argued that the Sogami-Ise theory 
contains a basic thermodynamic inconsistency, which when resolved leads to the usual 
DLVO potential. This, however, has failed to settle the issue and a number of papers 
are still being published on the thermodynamics of suspensions whose particles interact 
by the Sogami-Ise potential. 

A recently developed Digital Video Microscopy (DVM) has the potential of putting 
an end to this long standing debate. DVM provides the possibility of explicitly 
measuring the interaction potential between two macromolecules in suspension. For 
dilute bulk suspensions, the interaction potential between two spherical colloidal 
particles was found to be completely consistent with the DLVO theory ||119|| . A 
surprising result appeared, however, when a highly deionized suspension was confined 
between the glass plates [^, |7|, |176H173- It was discovered that when particles were 
close to a wall, the pair interaction potential developed a strong attractive component. 
The attraction was quite long ranged, comparable in its extent to the diameter of 
colloidal particle. 

It is well known that the glass-water interface has a significant negative charge due 
to dissociation of silanol groups. This charge density is comparable to that of colloidal 
particles. Furthermore, the absence of multivalent ions signifies that electrostatic 
correlations play only a marginal role. Under these conditions the Poisson-Boltzmann 
theory should work quite well. Indeed, the first numerical solution of the PB equation 
for two macromolecules inside a cylindrical pore has found an attraction between 
the two polyions [ p.78|| . The triumph of the PB theory was, however, short lived. 
Soon afterwards the mathematical proofs were published demonstrating that the PB 
equation is incapable of producing attraction between like-charged particles in a confined 
geometry ||179| , |18CI| , |181| , |182|| . The numerical calculation was fiawed, a consequence of 



the intrinsic difficulty of solving numerically non-linear partial differential equations with 
complicated boundary conditions. Indeed, it is extremely difficult to understand how 
a boundary can possibly change the interaction potential from repulsive to attractive. 
This is not to say that the presence of a boundary does not have a profound effect 
on the interactions. Consider one colloidal particle inside an electrolyte solution at 
fixed distance from a wall or an interface. In general, the interface is characterized by 
a dielectric discontinuity. For the moment, however, lets ignore this complication and 
suppose that the two sides have exactly the same dielectric constant. The interface, then, 
separates two half-spaces, one containing electrolyte and another electrolyte- free. The 
fact that all the charges are confined to one half-space strongly modifies the distribution 
of ions around the colloid. While the ionic cloud is spherically symmetric for colloidal 
particles far from the interface (in the bulk), it develops a strong asymmetry, which 
results in a net dipole moment of the colloid-electrolyte system. This dipole produces 
an electric field, 

E ~ ^ . (206) 
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which, because of the interface, can not be screened. If there are two macromolecules 
separated by a distance r along the interface, the electric field produced by one 
macromolecule interacts with the dipole moment p of the charge distribution induced 
by the second macromolecule. This leads to the effective interaction potential which is 
repulsive and falls off as 

w{r) ^ p ■ E ~ -3^ , (207) 

along the direction of an interface |p.83| , |184| , |185| , |186| , |187| , pl| . This is exactly the same 
kind of effective interaction found for 2d OCP, as seen in Section. |4.1| . 

It is clear that existence of an interface or a wall strongly modifies the interaction 
potential between two macromolecules. Instead of an exponential screening found in the 
bulk, the interactions along the interface fall-off algebraically as The potential, 

however, is still repulsive, and it is difficult to see how anything can modify this 
conclusion at the level of electrostatics. 

An interesting suggestion, which seems to account for the apparent attraction 
between the like-charged colloids near a wall, has been recently advanced by Squires 
and Brenner ||188|| . These authors attributed the attraction to non-equilibrium 
hydrodynamic flows which were not properly accounted for in experiment. 

While the hydrodynamics seems to be able to explain the apparent attraction 
between the colloids near a wall, it is not sufficient to explain the results of experiments 
in which colloids are sandwiched between two glass plates, since in this geometry the 



hydrodynamic attraction mediated by one wall is suppressed by the second wall |p.89 



Furthermore, hydrodynamics does not help to understand the long-lived metastable 
crystalline structures observed by Grier et al. | 177 , 19C] when a low density suspension 
is compressed against a glass plate. If the interactions between particles are effectively 
repulsive, once the constraint is removed the crystals should melt within seconds. 
Instead some crystalline regions are found to survive for as long as an hour. What 
is even more surprising is that the crystallites are actually three dimensional, extending 
far beyond the region were the pairwise surface-mediated attraction is found. This 
phenomenon is very similar to the voids observed by Ise et al.. 

The nature of confinement-induced attraction between the like-charged particles 
remains an open question. However, we must stress again that concentration on pair 
potentials when studying thermodynamic stability of colloidal suspensions is a serious 
oversimplification. The DLVO theory was proposed as an indicator of dynamical 
stability against flocculation, driven by short-ranged van der Waals forces. If the 
equilibrium structure of a colloidal suspension is in question, inter-colloid pair potential 
is not sufficient and the full free energy must be considered. We have already seen 
in Section ^3] that a large gain in solvation free energy obtained from the polyion- 
counterion interactions, strongly favors the phase separation of suspension into the 
coexisting high and low density phases. This tendency is opposed by the counterion 
condensation, which renormalizes the effective colloidal charge. The theory presented in 
Section suggest that colloidal suspensions should phase separate when C > 15.2. The 
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counterion condensation, however, prevents C from reaching this threshold. Nevertheless 
for highly charged colloidal particles, C can come very close to the critical value. 



Eq. (|125|) . This suggests that a deionized suspension of highly charged particles might 
actually be very close to criticality. This regime will be characterized by strong density 
fluctuations, which might appear as coexisting domains of voids and crystallites. 



8.2. Correlation-induced attraction 

DNA in aqueous solution is highly ionized due to dissociation of phosphate groups. This 
ionization results in one of the highest charge densities found in nature, one electron 
charge every 1.7 A. In spite of this huge charge concentration, over a meter DNA is 
packed into a nucleus of few micrometers. This efficient compaction is accomplished 
with the help of cationic proteins. The bacteriophages (viruses that infect bacteria) 
also use multivalent cations to package their DNA. Thus, the T7 bacteriophage head 
is 10~^ times smaller than the unpacked form of its DNA ||191|| . Furthermore, it is 
found that if the multivalent polyamines, known to exist in the host bacteria, are added 
to an in vitro solution containing DNA, the chains condense forming toroids very 
similar in size and shape to the ones found in vivo [H, 13. To prod uce condensation, 
multivalent counterions must somehow induce attraction between the different parts of 
[13, |l9l, |19|, |195, 



the DNA 



196 



The toroidal geometry is the result of the high 
intrinsic rigidity of the DNA molecule, which has the persistence length of = 500 A. 
If the compaction is done in such a way that the local radius of curvature, Tc, exceeds 
^p/27r the cost in elastic energy will be prohibitively high. Requirement that r^ > ^p/2tt, 
therefore, results in toroidal or spool-like condensates ||197| , |198|| . 

In eukaryotic cells the cytosol is traversed by a complicated network of 
microfilaments which are made of a protein called F-actin [^, |199|| . In spite of its high 
negative charge density F-actin, in the presence of multivalent counterions, agglomerates 
forming a network of bundles |pO0|| . Addition of monovalent salt screens the electrostatic 
interactions and re-dissolves the bundles lEOlll. What is the action of multivalent 



counterions which induces attraction between the like-charged macromolecules pH, E02 



203i |203| , 0, g |3|, 0, |T70|, |20|, ^O^, |208|, |0|, gig, |2n], |2T2|? To understand 



this we shall look at some very simple models p^ . 

Consider first two parallel polyions separated by a distance d inside a dilute solution 
containing a-valent ions. The polyions will be idealized as rigid lines of charge of length 
L = Zb. Each line has Z monomers of charge —q spaced uniformly along the chain. 
The solvent is a uniform medium of dielectric constant e. It is convenient to define the 
reduced polyion charge density as ^ = q^/eksTh [104]. A simple Manning argument 
then suggests that for ^ > 1 /a 

Z [ . 1 



1- , 



(208) 



a-ions condense onto each polyion. This is the lower bound on condensation, since the 
Manning argument does not take into account the correlations between the condensed 
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counterions. Nevertheless, even this simple estimate suggest that 88% of the DNA's 
charge should be neutralized by the divalent counterions. 

The associated counterions are free to move along the length of the DNA. We 
shall suppose that the only effect of condensation is the local renormalization of the 
monomeric charge from —q to (—1 + a)q. Lets define the occupation variables 0",^, with 
i = 1,2, Z and j = 1,2, in such a way that aij = 1, if a counterion is condensed at 
z'th monomer of the j'th polyion, and aij = otherwise. The occupation variables obey 
the constraint 

z z 

Y^aa = J2cTi2 = n, (209) 

i=l i=l 

where n is the number of condensed a-ions. The interaction energy between the two 
polyions is 

Z 2 

2e^^ r{i,3-i',3') ' ^ ^ 

1,1—1 ],f=i 

where the sum is restricted to (z, j) ^ and 



r{t,r, =h^\i- zf + (1 - 5,y)x2 (211) 

is the distance between the monomers located at and ii' ,]')■ Sjj' is the Kronecker 
delta, and x = d/h. The partition function is 

Q=Y^eM-m. (212) 

} 

where the prime indicates that the trace is done under the constraint of Eq. (|209|) . The 
force between the two polyions is 

This model is so simple that for polyions with not too high values of Z, the partition 
function can be solved explicitly For larger Z's the model can be easily simulated. 
In Fig. ^ we present the force as a function of separation for two polyions with Z = 20 
and n condensed divalent counterions. We see that in spite of the net like-charge, the two 
polyion-counterion complexes can attract each other at sufficiently small separations. 
Furthermore, we find that a critical number rii = Z/2a of condensed a-ions is necessary 
for the attraction to appear and that for monovalent counterions the interaction is always 
repulsive. 

What is the nature of this attraction? To understand this lets consider the limit 
of small separations between the polyions. The configurations which dominate the 
partition function for a; — > are the ones in which a condensed counterions on one 
polyion face the bare monomeric charges of the second polyion. In this limit the 
positions of condensed counterions on the two polyions become strongly correlated. 
The electrostatic energy of such configurations is 
E 2r^(l - ay , jZ - 2n)q^ 
ebx ebx 
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Figure 9. Force versus distance between polyions for Z = 20, a = 2, ^ = 2.283 
(corresponding to polymcthacrylatc) and n ~ 5, . . . , 10 (from top to bottom) in the 
Monte Carlo simulation [Q. Positive force signifies repulsion between the complexes, 
while the negative force implies existence of attraction. 



where n is the number of condensed counterions. The first term of Eq. (|214j) is the 
electrostatic energy due to attractive interaction between n condensed counterions 
and the bare monomers, while the second term is the repulsive interaction between 
Z — 2n uncompensated monomeric charges. We see that ij there are n > Z/2a 
condensed counterions the first term of Eq. ( pi4| ) dominates and the force becomes 
attractive for sufficiently short separations. For n < Z/2a the force is repulsive for all 
temperatures. This is the general mechanism of attraction between the two polyion- 
counterion complexes at any finite temperature. To minimize the electrostatic free 
energy, the positions of condensed counterions on the two polyions become correlated. 
For sufficiently short separations and number of condensed counterions exceeding 
the critical threshold ni, the correlation-induced attraction dominates the monopolar 
repulsion between the two complexes. 

The strength of the counterion correlations increase with a decrease in temperature. 
The state of maximum correlation is at T = 0. This, however, does not imply that the 
attraction between two polyions is maximum at zero temperature ||213|] . The reason for 
this is that the ground state configuration, which corresponds to the lowest electrostatic 
energy, is not, in general, the configuration which maximizes the attractive force. For 
low temperatures and short separations there are configurations which have force more 
attractive than the force in the ground state. For example, consider a polyion with 
Z = 6 and n = 3 condensed divalent counterions. The ground state corresponds 
to a staggered arrangement of counterions on the two polyions. The counterions on 
one polyion form the pattern H 1 1 — while on the second polyion they form a 
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Z^^^ = Z-anu = ^ ^ Z . (216) 



complimentary pattern — I 1 h- This leads to attraction between the complexes. 

However, it is easy to see that the configuration in which the counterions on the two 

polyions form the patterns + + H and h ++, has larger electrostatic 

energy, but leads to stronger attraction. At finite temperature the total force, being a 
weighted mean of forces associated with all the configurations can, therefore, become 
more attractive than the force at T = 0. This curious behavior, however, is confined 
to very small separations between the polyions. At larger distances the modulus of the 
attractive force is a monotonically decreasing function of temperature. 
If the number of condensed counterions exceeds 

= , (215) 

the interaction between the two lines of charge becomes, once again, repulsive. The 
riu corresponds to an overcharged configuration in which each polyion has the effective 
charge 

(a^l) 
2 

It is curious that exactly this kind of reentrant behavior is observed for the DNA 
condensates ||214| , p^l5| , |216| , |217| ]. In the absence of multivalent counterions, DNA inside 
solution has extended configuration. When the concentration of multivalent salt is 
slowly raised, there comes a point at which the gain in electrostatic energy due to 
polyion-a-ion association overcomes the entropic loss due to a-ion confinement in the 
vicinity of the polyions. The condensed monovalent counterions are then released into 
solution and are replaced by the polyvalent a-ions. After the critical number of a- 
ions is associated to the polyion, the interaction between the separate segments of the 
DNA becomes attractive. This plus the intrinsic rigidity of the molecule drives its 
condensation into toroidal bundles. As the concentration of multivalent ions in solution 
increases further, the DNA-a-ion complex becomes overcharged. When the effective 
charge of the polyion reaches Z^j^jr, the interaction between the segments of the DNA 
becomes, once again, repulsive and the bundles re-dissociate. Clearly this simple model 
is incapable of accounting for all the intricacies of the DNA condensation, nevertheless 
it sheds a lot of light on the role of electrostatic correlations in this interesting and 
important phenomenon. 

The calculations above were presented for a very idealized model of interacting 
lines of charge. It is quite simple to modify the theory to account for finite polyion 
diameter. This modification, however, does not significantly affect the predictions of 
the theory. Attraction appears at small separations between the polyion surfaces — 
about 7 A — after the critical number of a- ions is condensed onto the polyions [pl8|] . 
We find that for macromolecules of finite diameter less counterions are needed to induce 
attraction than for the two lines of equivalent charge density ||218|| . Furthermore, the 
charge-charge correlations along the polyion are of very short range [pl8| , |219|| , showing 
absence of any long-range order between the condensed counterions, contrary to the 
earlier speculations ||202| , |203| , 
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Figure 10. Two rodlike polyions of Z = 7. 



A very interesting experimental feature of rigid polyelectrolyte solutions containing 
multivalent counterions is that the correlation-induced attraction does not lead to 
phase separation [ 220|| . Instead rodlike polyions associate into bundles of well defined 
thickness. The precise bundle morphology depends on the persistence length of the 
polyions, all the bundles, however, tend to have a well defined cross-sectional diameter. 
What can account for this curious phenomenon? Correlation induced attraction should 



favor an unconstrained growth of bundles ||3^. In this respect the situation is very similar 
to that of an ionic crystal whose electrostatic energy is negative and is unbounded from 
bellow. What then cuts off the bundle size? The question is still not fully settled, but the 



indications are that the size of a bundle is controlled by the kinetics of its formation [p21 



To understand this better, lets, once again, consider our simple model of two lines of 
charge with ni < n < Z/a condensed a-ions pl3|| . The polyions are located on two 
parallel planes separated by a distance d. They are free to rotate in their respective 
planes. The relative angle between the two lines is 6*, with = corresponding to two 



parallel polyions. Fig. |10|. It is convenient to define the adimensional free energy as 

^=-^lng, (217) 
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Figure 11. Reduced free energy as a function of the angle 9 for Z = 9, rt = 4, 
a = 2, ^ = 2. The curves shown are for: x — 2 (full line), x — 1.3676 (dotted line), 
and X — 1.3 (dashed line). Notice that the free energies for 6 = and 9 = tt/2 are 
equal in the second case ^13| . 



where Q is the partition function for the two polyions at fixed separation and angular 
orientation and ^ is the Manning parameter. In Fig. |Tl| the reduced free energy is plotted 
as a function of 6 for various separations between the polyions. The free energy has 
two extrema at 6 = and 6 = 7t/2, corresponding to the parallel and the perpendicular 
orientation between the polyions. In general, in situations when the attraction appears, 
the perpendicular configuration has the lowest free energy at large separations between 
the polyions, while at small distances the parallel configuration minimizes the free 
energy of the model. However, in some cases, a reentrant behavior is observed. The 
perpendicular configuration is the global minimum for both large and small distances, 
while the parallel configuration minimizes the free energy at intermediate distances. 

At large separations the like-charged polyion-counterion complexes repel and are 
perpendicular to one another. As the distance between the polyions is decreases, the 
global minimum passes from 6 = tt /2 to 6 = 0, with 9 = becoming metastable. 
In order for polyions to form a bundle they must align, which means they have to 
overcome an activation barrier which separates the metastable minimum at 6 = 7r/2 
from the global minimum at 6 = 0. It has been suggested that the height of the barrier 
grows with the number of polyions which are already in the bundle p21|| . There comes 
a point when the barrier becomes so high that the thermal fluctuations are unable to 
overcome it and no new polyion can join the bundle. This puts an end to the bundle 
growth. 
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Figure 12. Two polyions with condensed counterions. Note that the net charge Zeff 
of one complex polarizes the counterions of the other complex. 



8.3. Counterion polarization 

The correlation- induced attraction between the polyions discussed in the previous 
section is short ranged, decaying exponentially with the separation between the polyions. 
The line of charge model is, however, too simple to account for the counterion 
polarization [p22|l . Consider for example a rigid cylinder-like polyion or a spherical 
colloidal particle with a layer of condensed counterions. Suppose that there is no 
overcharging so that only some fraction of the polyion charge is neutralized by the 
condensed counterions. The electric field produced by one complex polarizes the 
counterions of the other complex. Fig. |12[ This mechanism certainly provides an 
attractive component to the interaction potential |P23|| . However, can this induced 
attraction overcome the monopolar repulsion from the uncompensated colloidal charge 
Zeff7 To understand this lets consider two colloidal particles or globular proteins 
separated by a distance R inside an electrolyte solution with Debye length ^d- Each 
polyion has radius a, charge Z, and n condensed a-ions. The leading order interaction 
between the two macromolecules at separation i? > is of DLVO form 

VoLVoiR) = Zl^q^9\Kaf-^, (218) 

where Z^jj = Z — an The major corrections come from two effects. First, presence 
of polyions produces holes in the ionic atmosphere. The charge of a hole is minus 
the charge of the counterions excluded from the region of space occupied by the 
macromolecule ||117| , |118|| . For example, the charge of the hole produced by the polyion 



V2 in the ionic atmosphere polarized by the polyion Vi is 

« J-^^^ , (219) 

where the charge density at position R is given by the famihar Debye-Hiickel expression, 
P,iR) = -'-^^(R) = Zesq^^M^^a) ^ . (220) 



Use of the linearized Debye-Hiickel theory is justified by the fact that in Eq. ( p2(]| ) 
enters the renormalized effective charge. The electrostatic energy of the polyion-hole 
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interaction is obtained using the charging process in which q ^ Xq and 

rl 1 -2kR 

VhiR) = QhAR) J dM{R- \q) = -H'a^Zl^qH\Ka) , (221) 

at fixed K. The second major correction to the DLVO potential is the result of 
polarization of the condensed layer of counterions surrounding the polyion. Thus, the 
electric field produced by the complex V\ induces a dipole moment in the complex V^-, 
and vice-versa. The local surface charge density of V2 is 

a(r) = -(7_ + (T+e-^°'"^('") , (222) 

where 

and (/)(r) is the local electrostatic potential felt by the condensed counterions at position 



r. In the weak coupling limit F < 1 the exponential in Eq. ( |222| ) can be linearized 
yielding the surface charge distribution 

.(r) = Aa - ^ , (225) 
where Aa = o"+ — (T_ and the Gouy-Chapman length is 

'^'^ l-no^fiaq ^ ^ 

It is a straightforward calculation in electrostatics to show that the induced dipole 
moment of the complex V2 is 

2a^eE , 

The electric field produced by V\ at distance R from its center is 

E(R) = Z,^qQ{Kd) ^""''(1+^^-^)^ . (228) 

The electrostatic energy of the induced dipole of V2 inside the electric field produced 
by V\ is obtained using the charging process 

Vp = \ d\p{\q) ■ E(Ag) , (229) 
Jo 

where it is important to remember that Xgc{^q) = ^'^^gc{(1)- The same argument 
can be equally well applied to the interaction of V2 with the hole produced by Vi and 
to polarization of Vi by the electric field produced by 7^2- Summing up all of these 
contributions leads to 

W{R) = Zl^q'e\Kaf^ - Zl^q^ a' O^ko) h [j-] (230) 
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where the scaling function h{x) is 



(231) 



The correction to the DLVO potential is repulsive (hole dominated) for a/ \gc < 1-716... 
and is attractive (dipole dominated) for a/Xcc > 1.716.... We note, however, that it 
is always "doubly" screened and is, therefore, smaller than the leading order DLVO 
term for separations R > ||224|| . At shorter distances the approximations used to 
arrive at Eq. (|230|) fail, and the correlations between the condensed counterions must 
be explicitly taken into account [223]. For multivalent counterions, the electrostatic 
correlations result in short-ranged attraction, similar to the one discussed in the previous 
section. 



9. Conclusions 



We have explored the role of electrostatic correlations in systems ranging from classical 
plasmas to molecular biology. We saw how the positional correlations between the 
ions of an electrolyte can result in a thermodynamic instability. We also saw how the 
strong correlations between the polyions and the counterions lead to colloidal charge 
renormalization which stabilizes deionized suspensions against a phase separation. For 
two-dimensional plasmas electrostatic correlations are responsible for the metal-insulator 
transition. The critical behavior of superfluid ^He films, the roughening transition of 
crystal interfaces ||225|| , the melting of two-dimensional solids, and the criticality of XY 
magnets are all governed by the "electrostatic" interactions between the topological 
defects (charges). Nature has learned to take the full advantage of electrostatic 
correlations to efficiently package huge amounts of genetic material into tiny regions 
of space. 

Throughout this review we have come to rely on some simple models in order 
to understand complex physical phenomena. While these models are often sufficient 
to grasp the underlying physics, it is quite easy to push the models too far. 
This is particularly the case when one deals with specific structural properties of 
biomolecules ||226| , |227| , p28|| . As soon as the length scales on the order of few angstroms 
become important, approximation of water as a uniform dielectric medium is no longer 
sufficient \\229[ |230|| . Under these conditions reliance on simple models, which treat 
macromolecules and solvent as dielectrics, is probably no more than a wishful thinking. 
A careful path must be threaded between the simplification and the over-simplification. 
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